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ABSTRACT 

M17 is one of the youngest and most massive nearby star-formation regions in the Galaxy. It 
features a bright H II region erupting as a blister from the side of a giant molecular cloud (GMC). 
Combining photometry from the Spitzer Galactic Legacy Infrared Mid-Plane Survey Extraordinaire 
(GLIMPSE) with complementary infrared (IR) surveys, we identify candidate young stellar objects 
(YSOs) throughout a 1.5° x 1° field that includes the M17 complex. The long sightline through the 
Galaxy behind M17 creates significant contamination in our YSO sample from unassociated sources 
with similar IR colors. Removing contaminants, we produce a highly-reliable catalog of 96 candidate 
YSOs with a high probability of association with the M17 complex. We fit model spectral energy 
distributions to these sources and constrain their physical properties. Extrapolating the mass function 
of 62 intermediate-mass YSOs {M^, > 3 Mq), we estimate that >1000 stars are in the process of forming 
in the extended outer regions of M17. The remaining 34 candidate YSOs are found in a 0.17 deg^ 
field containing the well-studied M17 H II region and photodissociation region, where bright diffuse 
mid-IR emission drastically reduces the sensitivity of the GLIMPSE point-source detections. 

By inspecting IR survey images from IRAS and GLIMPSE, we find that M17 lies on the rim of 
a large shell structure '^0.5° in diameter (~20 pc at 2.1 kpc). We present maps of ^^CO and -'^^CO 
( J = 2 — > 1) emission observed with the Heinrich Hertz Telescope. The CO emission shows that the 
shell is a coherent, kinematic structure associated with M17, centered at u = 19 km s^^. The shell 
is an extended bubble outlining the photodissociation region of a faint, diffuse H II region several 
Myr old. We identify a group of candidate ionizing stars within the bubble. YSOs in our catalog 
are concentrated around the bubble rim, providing evidence that massive star formation has been 
triggered by the expansion of the bubble. The formation of the massive cluster ionizing the M17 H II 
region itself may have been similarly triggered. We conclude that the star formation history in the 
extended environment of M17 has been punctuated by successive waves of massive star formation 
propagating through a GMC complex. 

Subject headings: stars: formation — infrared: stars — infrared: ISM — H II regions — radio lines: 
ISM — dust, extinction 



1. INTRODUCTION 

Most star formation in the local universe is observed to 
occur in dense, massive clusters of hundreds or thousands 
of stars that are created by the gravitational collapse and 
fragmentation of massive cores within giant molecular 
clouds (GMCs). In the Galaxy, one of the best nearby 
laboratories for the study of star formation in the envi- 
ronment of a rich, massive cluster is M17, the Omega 
Nebula. Most of the observational attention given to 
M17 has focused on the bright "blister" H II region, ^5' 
in diameter, erupting from the side of the M17 molec- 

^ Department of Astronomy, University of Wisconsin at 
Madison, 475 N. Charter Street, Madison, WI 53706; 
povich@astro.wisc.edu 

^ Steward Observatory, University of Arizona, Tucson, AZ 85721 

^ Korea Astronomy and Space Science Institute, Chungnam Na- 
tional University, Daejeon, Korea 

University of Colorado, Space Science Institute, 1540 30th 
Street, Suite 23, Boulder, CO 80303-1012 

^ National Radio Astronomy Observatory, Charlottesville, VA 
22903 

^ Radio Astronomy Laboratory, University of California, Berke- 
ley CA 94720 

''' Department of Astronomy, University of Virginia, P.O. Box 
3818, Charlottesville, VA 22903-0818 

Spitzer Fellow, Harvard-Smithsonian Center for Astrophysics, 
60 Garden Street, Cambridge, MA 02138 



ular cloud, a GMC at a velocity of 20 km s ^ with a 
mass >3 x 10^ Mq (jLada II1976( ). The H II region di- 
vides the M17 molecular cloud into two components, and 
following the nomenclature of Wilson, Howe. & Balog3 
(jl99S) we will refer to the more massive, southern com- 
ponent as "M17 South" and the northern component 
as "M17 North" (see Fig. [1]). M17 has served both 
as an infrared (IR) spectral template for photodissoci- 
ation regions (PDRs) surrounding bright H II regions 
(jCesarskv et al.l 119961 : IVerstraete et al.l [19961) and as a 



prototype for H II region structure (Felli, Churchwell, 
& Massi 1984; Brogan & Troland 2001; Pellegrini et al. 
2007; Povich et al. 2007, hereafter PSC07). 

The ionization of the M17 H II region is provided 
by the rich cluster NGC 6618, which is deeply embed- 
ded in the gas and dust of the M17 cloud. Optical 
and near-IR photometric and spectroscopic studies have 
uncovered at least 16 O stars and over 100 B stars in 
NGC 6618 (Ogura & Ishida 1978; Chini, Elsasser, & 
Ncckel 1980; Lada et al. 1991; Hanson, Howarth, & 
Conti (1997); Hoffmeister et al. 2008, hereafter H08). 
Several water masers have been observed around M17 
(jCesarskv et al.|[T978D . and a prominent hypercompact 
H II region, M17 UC-1, lies on the interface of the H II 
region with M17 South (jFelli. Churchwell. fc Massil[l98l 
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ISewilo et al.ll2004) . Ml 7 has been searc hed repeatedly 
for ina.ssive stars with circ umstellar disks (jNielbock et alJ 
1200 It IChini et al.l l2004bl ). and thousands of IR-excess 
sources have been reported in and around the H II region 
(Lada et al. 1991; Jiang et al. 2002, hereafter J02; H08) 
High-resolution observations using the Chandra X-ray 
Observatory have provided a census of the young stellar 
population of M17 (Broos et al. 2007; hereafter BFT07). 
BFT07 found IR counterparts for 771 of the 886 X-ray 
sources in their sample, and they report that only ^^10% 
of these sources exhibit IR emission in excess of a stellar 
photosphere due to the presence of dust in circumstellar 
disks and/or infalling envelopes. Taken together, the ob- 
servations suggest that, while some fraction of the NGC 
6618 stellar population has ceased accreting, massive star 
formation is ongoing in and around the M17 H II region. 
Many authors have suggested that the birth of a second 
generation of massive stars has been triggered on the pe- 
riphery of the H II region by the ex pansion of the io nized 
gas into the M17 m olecular cloud (Cesar skv et al.iri978l : 
iFelh. Churchwell. fc Massi 1984 ; J02; HO^^^ 

The ex tended environment o f M17 has often been over- 
looked. lElmegreen fc Ladal ()1976l ) presented CO ob- 
servations of a large molecular cloud complex associ- 
ated with M17 at ?; = 20 km s~^ extending more 
than 1° to the Southwest, parallel to the Galactic mid- 
plane (Fig. [1]). They suggested that the entire com- 
plex is undergoing sequential star formation, begin- 
ning with the OB star s of N GC 6618 at the north- 
east tip. iJaffe fc Faziol (|1982l ) cast doubt on this se- 
quential star formation scenario by presenting 13 far- 
infrared sources within the extended complex that did 
not show a progression of decreasing age with in- 
creasing distance from Ml 7. Several high-resolution 
Galactic plane radio surveys in molecular or atomic 
line emis sion ha ve incl uded M17 ([ S anders et a l. 198ji; 
iMizimo fc Fukuil I2OOI iMcChire-Griffiths et all l2005h . 
The Massachusetts-Stony Brook Galactic Plane CO Sur- 
vey (| Sanders et al.l [l986l ). which provided {b,v) maps 
with 45" resolutio n that include M17 and its vicinity. 
The lSanders eTall ^1986) CO J = 1-0 maps clearly show 
bo th the M17 molecula r cloud and the complex identified 
by lElmegreen fc Ladal ()1976l ). and they reveal an addi- 
tional molecular cloud at ?; = 20 km s~^, centered at 
{l,b) = (15.9°, -0.7°). We name this second cloud MC 
G15.9-0.7 (Fig.[l]). 

Unbiased surveys of the Galactic plane reveal large 
structures that are missed by pointed observations with 
limited fields-of-view. Visual inspection of images from 
the Spitzer Space Telescope Galactic Infrared Mid-Plane 
Survey Extraordinaire (GLIMPSE) reveals a circular 
dust shell ^^0.5° in diameter outlined by faint, diffuse 
8 /im emission. M17 marks the Southwest rim of the 
shell (see Fig. [Ij. This structure is also apparent in 
Infrared Astronomy Satellite (IRAS) images at 60 /xm 
and 100 /im. The dust shell appears to extend from 
M17 to MC G15.9-0.7. Diffuse 8 /im emission is usually 
attributed to polycyclic aromatic hydrocarbons (PAHs) 
that fluoresce when illuminated by ultraviolet (UV) ra- 
diation. The shell morphology is simi lar to that of the 
IR bu bbles identified in GLIMPSE bv lChurchwell et all 
(|2006( 1. Such bubb les typica lly outhnc the PDRs of H II 
regions (|PSC07 : Watson et al., 20 08). The dust shell near 
M17 was not selected for the GLIMPSE bubbles cata- 



log because it is partially obscured by the bright diffuse 
emission from M17 and lacks a well-defined inner rim. 
Nevertheless, we will demonstrate that the shell is an ex- 
tended bubble outlining the PDR of a large, faint H II 
region associated with M17. We call this structure the 
M17 extended bubble (M17 EB; Fig.[l]). 

While the screen of cold, dense molecular gas in M17 
South probably prevents the stellar winds and ionizing 
radiation from the OB stars in NGC 6618 from influenc- 
ing the extended molecular cloud complex to the South- 
west, on the Northeast side of M17 there is evidence that 
M17 EB is interacting with both MC 15.9-0.7 and M17 
itself. The M17 H II region surrounds a central cavity 
that is filled by hot, diffuse plas ma observed by C/ian- 
dra (jTownslev et al.ll2003HPSC07f ). i?05'^r observations 
covering a larger area show soft X-ray emission spilling 
out of the H II region and apparently filling part of M17 
EB; this X-ra y emission is spatially coincident with faint 
Ha emission (|Dunne et a l. 2003). A search for SIMBAD 
objects in the interior of M17 EB reveals several stars 
with OB spectral types. We will show that these stars 
represent the most massive members of the cluster ion- 
izing M17 EB and that this cluster is both distinct from 
and older than NGC 6618. We refer to this progenitor 
cluster as NGC 6618PG. 

The distance to Ml 7 has been somewhat disputed in 
recent years. Most estimates have been based upon 
optical/near-IR photometry of th e most luminous stars. 
iHanson. Howarth. fc Contil ()1997l ) reported a distance of 
1-3-0 2 kpc, a significant departure from the kinematic 
distance of 2.3 kpc. iNielbock et al.l ()2001l) subsequently 
empl oyed a revised extinction law (jChini fc Wargaul 
Il998f ) and derived a distance of 1.6 ± 0.3 kpc. The value 
of 1.6 kpc has since been widely adopted. PSC07 bal- 
anced the flux in the integrated IR SED of Ml 7 with 
the bolometric luminosity of all known ionizing stars in 
NGC 6618 and reported a distance of 1.6^0^ kpc. This 
luminosity distance assumed both that all of the ma- 
jor ionizing stars in the region were known and, like the 
previous spectrophotometric distance estimates, that the 
majority of the known ionizing stars were not in unre- 
solved binary or otherwise multiple systems. Recently, 
|H08| showed that the most luminous stellar system in 
M17, the 04-1-04 binar y star called Kleinmann's Anony - 
mous star or GEN 1 (jChini. Elsasser. fc Neckel Il98tf). 
is itself composed of two spectroscopic binary systems 
of nearly equal-mass components. This discovery, along 
with new distan ce moduli for other massive stars in NGC 
6618, lead H08 to revise the M17 distance once again, to 
2.1 ± 0.2 kpc, in agreement with the kinematic distance. 
In this work, we will generally adopt the 2.1 kpc distance, 
but we will also consider 1.6 kpc where the distance un- 
certainty could significantly impact our results. 

We take advantage of the large area surrounding M17 
observed by the Spitzer Galactic plane surveys to in- 
vestigate star formation in the extended environment of 
M17, emphasizing the newly-discovered structures to the 
Northeast. We have combined GLIMPSE with longer- 
wavelength IR survey data from both the Multiband 
Imaging Photometer for Spitzer (MIPS) Galactic plane 
survey (MIPSGAL) and the Midcourse Space Experiment 
(MSX). We use these datasets to investigate the popula- 
tion of young stellar objects (YSOs), identified by their 
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Fig. 1. — GLIMPSE full-color image of the M17 field {red: 8.0 nm, orange: 5.8 fim, green: 4.5 fxm, and blue: 3.6 pm). 
A schematic diagram showing the approximate boundaries of large-scale structures men tioned in the text is overlaid. The 
molecular clouds shown in white and gray have been identified from the CO survey data of lSanders eraD ^M)- 



IR excess emission, over a large (~ 1.5° x 1°) area sur- 
rounding M17. Our investigation includes the full spatial 
extent of the M17 molecular cloud, M17 EB, and MC 
G15.9-0.7. We have followed up the infrared survey data 
with high-resolution ^^CO and ^'^CO J — 2-1 observa- 
tions. In §2 we summarize our observations and data 
reduction. We describe our methods for selecting and 
characterizing YSO candidates in §3, and we discuss the 
extended young stellar population of M17 in relation to 
the molecular cloud structures in §4. In §5, we synthesize 
these results into a picture of propagating massive star 
formation in M17, and we summarize our conclusions in 
§6. 

2. OBSERVATIONS AND DATA PROCESSING 
2.1. IR Galactic Plane Survey Data 
2.1.1. GLIMPSE 

As part of the Spitzer Lega cy Science Program, 
GLIMPSE (|Beniamin et al.ll2003f) imaged the Galactic 
plane from \l\ = 0°-65°, |6| < 1° using the four mid- 
IR bands (3.6, 4.5, 5.8, and 8.0 /im) of the Infrared 
Array Camera (IR AC) on the Spitzer Space Telescope 
(jFazio et al.ll2004f ). High-resolution (1.2" pixels) image 
mosaics were created by the GLIMPSE pipeline^ from 
Basic Calibrated Data (BCD) image frames processed by 
the Spitzer Science Center (SSC). The GLIMPS E point- 
source extractor, a modified version of Daophot ((StetsonI 

^ Details of the data processing ca n be found at 
|http : // www .astro.wisc.edu/ glimpse /docs . html | 



was used to extract point sources from the individ- 
ual image frames. The 5a point-source detection limits 
of GLIMPSE are nominally 0.2, 0.2, 0.4, and 0.4 mJy 
for the [3.6], [4.5], [5.8], and [8.0] bands, respectively, 
and higher in regions of bright diffuse emission. Flux 
densities of sources detected at greater than 5a at least 
twice in one or more of the 4 IRAC bands were included 
in the GLIMPSE Point Source Archive. The highly- 
reliable GLIMPSE Point Source Catalog is a subset of the 
Archive. Both the Archive and the Catalog incorporated 
JHK s flux densities from the 2MASS point source cat- 
alog (iSkrutskie et al.l [2006( 1 to produce final sourcelists 
containing flux densities from 1-8 /xm in 7 near- to mid- 
IR bands. For this study, we use the GLIMPSE v2.0 
Point Source Archive and the GLIMPSE v3.0 image mo- 
saics. The majority of the Archive sources used in our 
analysis will also appear in the highly-reliable Catalog 
because of our additional requirement that each source 
must be detected in at least 4 bands (see §3.2). We use 
the Archive instead of the catalog because the Archive 
identifies sources with near-saturated fiuxes that can 
later be used as lower limits and contains more sources 
in crowded regions. Mosaics produced from 2MASS data 
provide near-IR images that complement GLIMPSE. 

2.1.2. MIPSGAL 

We downloaded 24 /im and 70 /im post-B CD mosaics 
of the MIPSGAL survey (iCarev et al.ll2009f l from SSC. 
The 24 /im mosaics have 2.4" pixels and the 70 /im mo- 
saics have 4.8" pixels. MIPSGAL point-source catalogs 
were not yet available at the time of this study, so we ex- 
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tracted 24 /xm point-source fluxes using the GLIMPSE 
extractor. Because the GLIMPSE pipehne software is 
not optimized for MIPSGAL data, this procedure gener- 
ates the foflowing types of false source detections: noise 
peaks, small resolved structures, and artifacts from the 
point-spread function wings of bright sources. As a first 
cut to reduce the number of false sources in our analy- 
sis, we removed all 24 /im sources with larger than 25% 
error in measured flux density [F/SF < 4) from the 
sourcelist. We then spatially correlated (bandmerged) 
the 24 fj,m sources with the GLIMPSE Archive sources. 
If an Archive source was located within 2.4" of the central 
position of a 24 /xm source, the source was considered to 
have a 24 /im detection and the [24] flux was appended to 
the Archive fluxes. If more than one Archive source was 
located within 2.4" of a 24 /im source, the [24] flux was 
appended as an upper limit for each correlated Archive 
source. The correlation procedure favors "true" 24 /xm 
point sources because these are more likely to have a 
GLIMPSE counterpart, but occasionally a false source or 
false correlation makes it into the bandmerged list. The 
flnal 8-band (combined 2MASS, GLIMPSE, and MIPS- 
GAL [24]) sourcehsts were inspected by eye to reject any 
remaining suspicious sources (see §3.2). 

2.1.3. Aperture Photometry 

We did not attempt to automatically extract MIPS- 
GAL [70] point source fluxes. Only ~10 obvious 70 /im 
unresolved sources are visible in the BCD mosaics within 
our analysis region, so we extracted the [70] fluxes of in- 
teresting sources manually, using aperture photometry. 
Visual inspection of the GLIMPSE mosaics revealed ^10 
compact red sources that are resolved by IRAC and hence 
not in the Archive. Such sources may be luminous YSOs 
that heat ambient dust in molecular clouds out to a few 
pc distance. One pc corresponds to 78 1.2" GLIMPSE 
mosaic pixels at 2.1 kpc. Our aperture photometry pro- 
cedure used the Funtools^° package with the SAOIm- 
age DS9^^ image display program to extract fluxes from 
apertures deflned by eye. For each source extracted, we 
used the same target apertures for the 7 2MASS-I-IRAC 
bands and increased the aperture size as needed for the 
2 MIPSGAL bands. We varied the position of the back- 
ground sampling a minimum of 3 times to estimate the 
contribution of the background subtraction to the flux 
uncertainty. Although the background subtraction gen- 
erally dominated the uncertainty on the extracted fluxes, 
our uncertainty estimates also included photon counting 
statistics (prior to the correction for Zodiacal emission) 
and the uncertainty introduced by applying the IRAC 
extended aperture corrections^^ to sources with effective 
radii > 9". 

2.1.4. MSX 

The Spirit III instrument on board the M SX satellite 
surve yed the Galactic plane in four IR bands (jPrice et al.l 
I200H ): A (isophotal wavelength 8.28 /xm), C (12.13 /fm), 
D (14.65 /xm), and E (21.3 /im). The 5a point-source 
detection limits for MSX are 100, 1100, 900, and 200 

See http://hea-www.harvard.edu/ RD/funtools/ ds9.htmll 

See http://hea-www.harvard.edu/RD/ds9/refi 

See http://ssc.spitzer.caltech.edu/irac/calib/extcal/ 1 



mJy for A, C, D, and E bands, respectively. The sen- 
sitivity of the MSX A and E bands is just below the 
saturation limits of IRAC and MIPS, so these bands can 
be used to replace the fluxes of point sources saturated 
at IRAC [8.0] and MIPS [24], respectively The spatial 
resolution of the MSX images is ~18.3", so confusion is 
an issue when correlating MSX Catalog^"^ sources with 
the GLIMPSE Archive. After selecting YSO candidates 
based on the combined GLIMPSE and MIPS [24] pho- 
tometry (see §3.1), we correlated only the brightest satu- 
rated and near-saturated YSO candidates with the MSX 
catalog, using a 6" correlation radius. We did not use 
the IRAS Point-Source Catalog for this work because its 
low resolution causes too much confusion with multiple 
GLIMPSE sources. 

2.2. CO Observations with the Heinrich Hertz 
Submillimeter Telescope 

We mapped the M17 region in the J = 2 — > 1 transi- 
tions of i^CO and ^^CO with the 10-m Heinrich Hertz 
Telescope (HHT) on Mt. Graham, Arizona between 2008 
February and 2008 June. The receiver was a dual polar- 
ization superconductor-insulator-superconductor (SIS) 
mixer system employing single-sideband (SSB) mixers 
with outputs for both upper and lower sidebands, each 
with a 4-6 GHz IF band. The ^^CO J = 2 ^ 1 line 
at 230.538 GHz was placed in the upper sideband and 
the i^CO J = 2 ^ 1 line at 220.399 GHz in the lower 
sideband, with a small offset in frequency to ensure that 
the two lines were adequately separated in the IF band, 
since the sideband rejection was typically 17-20 dB. The 
spectrometers, one for each of the two polarizations and 
the two sidebands, were filter banks with 256 channels 
of 1 MHz width and separation. At the observing fre- 
quencies, the spectral resolution was 1.3 km s~^ and the 
angular resolution of the telescope was 32" (FWHM). 
A set of 10' X 10' fields were mapped with on-the-fly 
(OTF) scanning in RA at 10" per sec with row spacing 
of 10" in declination, over a total of 60 rows. Each field 
required about 2 hours of elapsed time. System tem- 
peratures were calibra ted by the standard am bient tem- 
perature load method (jKutner fc Uhch|[l981[ ) after every 
other row of the map grid. Atmospheric conditions were 
generally clear and stable, and the system temperatures 
were usually nearly constant at Tgys ^ 200 K (SSB) for 
each 10' x 10' field. 

Data for each polarization and CO isotopomer were 
processed with the CLASS reduction package (from the 
University of Grenoble Astrophysics Group), by remov- 
ing a linear baseline and convolving the data to a square 
grid with 10" grid spacing (equal to approximately one- 
third the telescope beamwidth). The intensity scales for 
the two polarizations were determined from observations 
of W51D (W51-IRS2) made just before the OTF maps. 
The gridded spectral data cubes w ere processed with the 
Miriad software package (ISault e t al. 1995) for further 
analysis. The two polarizations were averaged for each 
sideband, yielding images with RMS noise per pixel and 
per velocity channel of 0.15 K-T^ for both the ^^CO and 
^■^CO transitions. 

We mapped a total of 26 10' x 10' fields, arranged in a 
semi-regular pattern to cover the main M17 H II region 

See |http://irsa.ipac. caltech.edu/applications/Gator/ 1 
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and M17 EB. These 26 fields were combined onto a com- 
mon spatial grid and converted to Galactic coordinates 
for comparison with the GLIMPSE images. 

3. YSO SELECTION AND CHARACTERIZATION 

An overview of the target field of the GLIMPSE survey 
searched for YSOs is presented in Figure [TJ with the ap- 
proximate boundaries of the M17 H II region, the bubble 
M17 EB, and t he m ajor molecular cloud structures from 
ISanders et al.l (1986) overlaid. The top of the image is 
the Galactic mid-planc (5 = 0°), and the jagged bottom 
edge is the boundary of the GLIMPSE survey between 
b = —1° and —1.2°. Most of this image is occupied by 
the molecular cloud complex at w = 20 km s~^ that in- 
cludes M17. A smaller cloud, MC G14.9-0.4 associated 
with compact, bright diffuse mid-IR emission, lies in the 
background at w = 60 km s~^. 

We identify candidate YSOs in this image based on 
their IR excess emission, with the goal of determining 
the YSO population associated with M17. The location 
of M17 near the Galactic mid-plane and only 15° from 
the Galactic center produces a long sighthne through the 
Galaxy in the mid-IR that passes through multiple spi- 
ral arms; hence we expect that the GLIMPSE Archive 
contains many YSOs or other objects with similar IR 
colors that lie at different distances, unassociated with 
M17. The average extinction produced by the molecu- 
lar clouds in the target field is too low to significantly 
reduce the number of background sources detected in 
the mid-IR (see §4.1.1). Because M17 is bounded in 
longitude by potentially star-forming molecular clouds 
and in latitude by the Galactic mid-planc and the sur- 
vey edge, there is no portion of the GLIMPSE image 
in Figure [T] that could be used as an off-source region 
to estimate the contamination in our sample from fore- 
ground or background sources (see §3.3). Fortunately, 
the region symmetric to M17, reflected above the Galac- 
tic mid-plan e, lacks obvious CO emission marking molec- 
ular clouds (jSanders et al.lll986HM"izuno k Fukuill2004f) 
and is correspondingly free of the bright diffuse 8.0 /xm 
emission in GLIMPSE that usually identifies H II re- 
gions. We selected a circular control field of 0.5° radius 
centered at {l,b) = (15.08°, 0.61°). Since this region has 
the same Galactic longitude and absolute latitude occu- 
pied by M17 but exhibits no obvious signs of massive 
star formation, it can be used to address the question of 
what the Galactic plane would look like if M17 were not 
present. We analyzed the control field in parallel with 
the M17 target field, using the same procedure to select 
candidate YSOs from both fields. 

3.1. Fitting YSO Model SEDs to Broadband IR Fluxes 

Our principal tools for identifying and characterizing 
candidate YSOs from this rich datas e t wer e the grid 
of YSO models from iRobitaille et all ()2006l . hereafter 
RW06) an d the spectral energy distribution (SED) fit- 
ting tool of lRobitaille et al.l ()2007l . hereafter RW07). The 
model grid consists of 20,0 00 2-dimensional Mon te Carlo 
radiation transfer models (jWhitnev et al.ll2003allbl . I2004f) 
spanning a complete range of stellar mass and evolution- 
ary stage and output at 10 viewing angles (inclinations), 
so the model fitting tool actually has 200,000 SEDs to 
choose from. The model fitting tool uses a fast x^" 
minimization linear regression algorithm (RW07). We 



can robustly distinguish between YSOs and reddened 
photospheres of main-sequence and giant stars because 
YSOs require a thermal emission component from cir- 
cumstellar dust to reproduce the shapes of their mid- 
IR excesses. The Monte Carlo radiative transfer models 
upon which the fitting tool relies have been tested ex- 
tensively by successfully fitting the SEDs of numerous 
well-characterized YSOs in Taurus (RW07) . The models 
and fitting tool have also been emplo yed to analyze the 
YSO populations of the Eagle Nebula (jlndebetouw et al] 
|2007() . several other Galactic massive st ar formation re- 
gions observed as part o f GLIMPSE (| Shepherd et al.l 
I2007t iWatson et alT 1200 8'). and the Large Magellanic 
Cloud (fWhitnev et al.l l2008'l. 

The radiation transfer technique employed in the 
RW06 models propagates "photons" from the central 
source through the circumstellar environment. The mod- 
els solve for the temperature structure of the circumstel- 
lar material and include absorption and re-emission by 
dust along with photons produced by dis k accretion and 
back war ming of the stellar photosphere (jWhitnev et al.l 
I2003allbl . I2004D . When we fit models to observed YSOs, 
we essentially use Monte Carlo radiation transfer to 
look through the veil of circumstellar dust, placing 
the central star on the H-R diagram by interpolating 
b et ween pre-main-sequence ( PMS) evolutionary tracks 
(IBernasconi fc Maedeii Il996t iSiess. Dufour. fc Forestinil 
I2000D . Like all complex models, the YSO models 
have built-in assumptions, degeneracies, and a multi- 
dimensional parameter space to explore. In general, fit- 
ting these models to real data results in many differ- 
ent sets of parameters that can describe a YSO almost 
equally well. Rather than attempt (futilely) to find a 
unique solution, we accept a range of well-fit models to 
investigate how well we can constrain the physical prop- 
erties of each YSO. The distance to the modeled sources 
places an important external constraint on the allowable 
luminosities of the well-fit YSO models. For M17, we 
used the distance range of 1.6-2.3 kpc, effectively build- 
ing the distance uncertainty into our constraints on the 
YSO properties from the model parameters. The number 
of acceptable models also decreases rapidly with increas- 
ing range of IR wavelengths from available data included 
in the SED. 

The steps of our process for selecting candidate YSOs 
are outlined in Table [1] along with the number and den- 
sity of sources in both the target and control fields at each 
step. We began by removing sources with SEDs that are 
consistent with stellar photospheres. We fit all sources 
that are detected in N^ata > 4 of the 8 near- through 
mid-IR bands in the GLIMPSE Archive combined with 
MIPSGAL 24 /im fluxes (126,385/352,225 sources in the 
target field and 71,391/191,010 sources in the control 
field; Ta ble □ with SEDs from 78 53 model stellar atmo- 
spheres (i Brott fc Hauschildtll2005f) included in the RW07 
fitting tool. Interstell ar extinction based upon the mid- 
IR extinction law of llndebetouw et ahl (j2005) was in- 
cluded in the model fits, so even highly reddened stars 
(we allow a maximum Ay of 30 mag when fitting stellar 
photosphere SEDs) returned good fits. We considered all 
sources with a best-fit x^, normalized by the number of 
flux datapoints used in the fit, satisfying /N^ata < 5 
to be well-flt by stellar photospheres. We determined 
this threshold value by visual inspection of the best-fit 
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TABLE 1 

Source Counts in the YSO Search 





M17 Target Field 


Control Field 




(1.67 deg2) 


(0.78 deg2) 


Sources 


Number 


Density (deg ) 


Number 


Density (deg ) 


In GLIMPSE Archive 


352,225 


211,000 


191,010 


245,000 


Fit with SED models (A^data > 4) 


126,385 


75,700 


71,391 


91,500 


Well-fit by stellar photosphere SEDs (x^/A^data < 5) 


124,768 


74,700 


70,628 


90,500 


Poorly-fit by stellar photosphere SEDs (x^/A^data > 5) 


1,617 


970 


763 


980 


Well-fit by YSO SEDs (xV^^data < 5) 


979 


590 


370 


470 


In final sample of YSO candidates 


406 


240 


106 


140 


In M17 YSO candidate subsample^/control subsample'' 


195 


205 


25 


100 



^ Irregularly shaped region with area 0.95 deg^ (Sec §4). 

^ Square box centered at (i, b) = (15.1°, 0.7°) with area 0.25 deg^. 



SEDs for a variety of sources. To avoid biasing the fits for 
sources where the flux uncertainties have been underesti- 
mated, before fitting any models we conservatively reset 
the uncertainties to a floor value of 6F/F = fO% for 
all sources in the GLIMPSE Archive with dF/F < 10%. 
For the MIPS [24] fluxes, we used an uncertainty floor of 
15%. 

We then took the 1,617 sources in the target field and 
763 sources in the control field that were poorly fit by 
stellar photospheres and fit them a second time using 
the RW06 YSO models. Again, we considered sources 
with a best-fit x^/^data < 5 to be well-fit, and these 
formed the basis of our YSO sample. The vast major- 
ity of sources that were well-fit neither by stellar photo- 
spheres nor YSO models are actually stars with saturated 
fluxes or poor photometry (often in the [24] band where 
our photometry is less robust) or signs of variability (a 
mismatch between 2MASS and GLIMPSE photometry 
due to the different epochs of the observations). Eleven 
sources in the target field appeared to be good YSO can- 
didates with marginal fits. Several of these had satu- 
rated fluxes in one or more bands, and one (G015.6653- 
00.4989) exhibited excess emission in the [4.5] band 
similar to previou sly-studied sources associated with 
molecular outflows (ISinith et al.ll2006l: [Pavis et al.ll2007l : 
iShepherd et~al] l2007t iCveanowski et all I2008D . These 
were moved to the "well-fit" sample after the question- 
able flux measurements were replaced with upper or 
lower limits. 

3.2. A Highly Reliable Sample of Candidate YSOs 

Our primary goal was to produce a highly reliable 
sample of candidates YSO. As shown in Table [1] only 
979/126,385 or 0.77% of sources fit with SED models in 
the M17 target field were fit wcU by YSO models but 
could not be fit well by stellar photospheres. This frac- 
tion was significantly lower for the control field, with 
370/71,391 or 0.5% of sources (the counting uncertainty 
is 0.02%). While the model fitting process performed the 
service of discarding the >99% of sources that we were 
not interested in studying, simply being well-fit by YSO 
model SEDs was a necessary but insufficient criterion for 
a source to be included in the final sample. As a last 
step, we culled sources from the sample by inspecting all 
of the SED fits by eye, in many cases returning to the 
GLIMPSE and MIPSGAL images to confirm visually the 
existence of a valid source in each band with a reported 



detection. Two common reasons to reject a well-fit YSO 
candidate necessitated this final inspection: 

1. The IR excess emission that prevented the source 
from being fit well by stellar photosphere models 
occurs only in the IRAC [8.0] band, with no detec- 
tion at MIPS [24]. Such an SED can be produced 
by poor source extraction in regions with highly 
structured diffuse emission or a Malmquist bias af- 
fecting faint sources in the GLIMPSE sourcelists. 
The Malmquist bias occurs because when a source 
is selected for the Archive on the basis of a 5(7 de- 
tection in any combination of bands, the fluxes of 
the remaining bands can be entered in the Archive 
even if they are at lower confidence levels. Stars 
are typically faint at 8 /im, and often a noise peak 
or a diffuse background feature can be extracted 
from the position of a star observed in the other 
bands, causing an artificially high measured [8.0] 
flux. While a few of the sources exhibiting this SED 
may actually be YSOs, we make the conservative 
decision to discard all questionable candidates. 

2. A suspicious IR excess is observed in the MIPS 
[24] band only. The difference in resolution be- 
tween MIPS and IRAC increases the chance of 
a false match when correlating sources between 
MIPSGAL and GLIMPSE. This often happens 
near bright, saturated 24 /im sources where the ex- 
tended wings of the MIPS point-spread function 
can be extracted as individual point sources and 
matched with an overlapping GLIMPSE source. 
This creates a spurious 24 /im excess attributed 
to the GLIMPSE source. 

Both of the above issues constitute important caveats 
to consider if searching for YSOs possessing "transition 
disks" with large inner holes that might have been cleared 
out by planets, since such sources exhibit IR excess emis- 
sion only at longer wavelengths. 

The brightest sources in our sample also appear in the 
MSX point source catalog, so we used available MSX 
fluxes to fill in the SED between IRAC and MIPS and 
replace saturated IRAC [8.0] and MIPS [24] fluxes, where 
possible. Bright sources with very red SEDs from K- 
[8.0] that are isolated from other IR-excess sources are 
candidate evolved stars on the asymptotic giant branch 
(AGB). Carbon stars, for e xample, have mean absolute 
K-heoid magnitudes of ^^6.8 (IWallerstein fc Knapd[l998[) 
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Fig. 2. — Example plots of best-fit model SEDs to available near- and mid-IR photometry. Two sources well-fit by AGB star 
spectral templates {left panels) are compared to two sources well-fit by YSO models {right panels) with similar brightnesses. 
Triangles denote datapoints with saturated fluxes that were used only as lower limits for the model fitting. The best-fit SED to 
the KW Object is plotted in the upper-right panel. 



and are often en shrouded by dusty env elopes generated 
by strong winds (|Guandahm et al.i2006f) : such stars have 
similar SEDs to YSOs from K to [8.0] and, even when 
located on the far side of the Galaxy, will be among 
the brightest sources observed in GLIMPSE. The RW07 
fitting tool incorporates 75 AGB star SED templates 
derived f rom Infrared Space Observatory {ISO) spec- 
trosc opy (jSvlvester et al.llT999l : iHonv. Waters, fc TielensI 
|2002|) . After fitting bright sources with the AGB star 
templates, we removed 30 and 15 likely AGB stars from 
the target and control samples, respectively. 

Example SEDs of two probable AGB stars and two 
bright candidate YSOs from the target field are plot- 
ted in Figure [21 The most distinctive features of the 
AGB star SEDs are the precipitously red if-[8.0] spec- 
tral index and the turnover to a Rayleigh- Jeans spectrum 
at wavelengths longward of ~24 ^m. The YSO models 
that fit the SEDs of these bright sources tend to be mas- 
sive stars with optically thick disks but without substan- 
tial infalling envelopes. Such objects are relatively rare, 
the best exam ple in M17 being the Klei nmann- Wright 
(KW) Object (iKleinmann fc Wrighi3ll973) . a suspected 
Herbig Be star (jChini et al.l l2004a[) . The SED of the 
KW Object is plotted in the upper-right panel of Fig- 
ure [2l Several candidate YSOs, including the one shown 
in the lower-right panel of Figure [H are detected at 70 
fj,m by MIPSGAL, but no candidate AGB stars are ob- 
served at 70 ^m. Generally, YSO disks and envelopes 
contain cool (T^ ~ 30 K) dust that becomes optically 
thin at relatively large radii from the central star, while 
the dust distribution in AGB stellar winds has a much 



smaller effective radius, giving higher effective tempera- 
tures {Td > 200 K). It is therefore reasonable to expect 
that the SEDs of AGB stars will peak at shorter wave- 
lengths than those of YSOs, but this difference cannot 
be observed without data at wavelengths of 24 /im or 
longer, and there may be exceptions to this general rule. 

Visual inspection of the M17 target field revealed 6 
bright, compact resolved IRAC sources that have SEDs 
well-fit by YSO models. These sources are noted in 
Tables EH (see §4.1.2). Four of the 6 are MSX point 
sources, and their SEDs are suggestive of a PAH spec- 
trum. The model fits to these sources were improved 
dramatically by assuming PAHs contributed significantly 
to the emission in the 4 bands with the strongest PAH 
features, IRAC [5.8] fc [8.0] and MSX A & C. The RW06 
YSO models did not include emission from PAHs; hence 
to fit sources with potentially PAH-dominated spectra we 
set these 4 bands as upper limits. Using the fluxes of the 
well-fit models to estimate the continuum levels in the 
bands affected by PAHs, we find that PAHs contribute 
anywhere from 30% to 90% of the broadband fluxes. The 
YSO models for 3 of these 4 sources indicate central stars 
with photospheric temperatures of Tgs 2 x 10^ K, hot 
enough to produce the UV radiation necessary to excite 
PAH transitions (PSC07, and references therein). The 
fourth source, F53 in Table [31 was fit by few high- Toff 
models. This object appears to be associated with the 
molecular cloud MC G14.9-0.4 (see Fig. [J), giving it a 
near kinematic distance of 4.8 kpc at the cloud velocity 
of 60 km s~^. This is far beyond the maximum distance 
of 2.3 kpc given to the model fitting tool, and at the 
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TABLE 2 

Candidate YSOs in the Extended Environment of M17 



Index 


I 


h 


[4.5] 


[8.0] 


[24] 




(M©) 


log(L' 


tot/Lq) 


Evolutionary 


No. 


(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


Min 


Max 


Min 


Max 


Stage^ 


The Bridge 


El 


14.7207 


-0.4859 


6.4 


6.0 


3.3 


6 


6 


3.1 


3.2 


II 


E2 


14.7392 


-0.4736 


12.0 


10.4 


5.7 


0.3 


4 


0.1 


2.2 


Amb. 


E3 


14.7762 


-0.4657 


7.4 


6.8 


3.9 


6 


6 


3.1 


3.2 


II 


E4 


14.7776 


-0.4900 


9.7 


9.5 


5.8 


1 


4 


1.3 


2.1 


Amb. 


E5 


14.7963 


-0.4565 


12.3 


11.2 


6.9 


1 


4 


1.3 


2.5 


II 


E6 


14.8039 


-0.5135 


7.4 


7.0 


4.3 


3 


6 


2.1 


2.9 


II 


E7 


14.8214 


-0.5136 


7.5 


6.5 


2.8 


4 


8 


2.4 


3.4 


II 


E8 


14.8217 


-0.4782 


11.3 


11.1 


6.8 


0.6 


4 


0.6 


1.7 


O/I 


E9 


14.8291 


-0.5279 


9.2 


8.4 


4.8 


3 


5 


1.8 


2.6 


II 


ElO 


14.8462 


-0.4890 


12.7 


11.9 


7.0 


0.2 


4 


0.1 


1.7 


Amb. 


The Outer Regions of the M17 CMC 


Ell 


14.8338 


-0.6974 


9.6 


8.8 


6.0 


3 


5 


1.7 


2.6 


II 


E12 


14.9743 


-0.5079 


11.7 


9.8 


> 3.3 


0.3 


4 


0.4 


2.3 


Amb. 


E13 


15.1289 


-0.4958 


11.5 


9.0 


3.2 


0.2 


5 


1.0 


2.7 


O/I 


E14 


15.1426 


-0.4883 


11.4 


9.8 


> 4.1 


0.2 


5 


0.5 


2.8 


II 


E15 


15.1473 


-0.5299 


9.2 


8.3 


3.5 


1 


5 


1.5 


2.5 


Amb. 


E16 


15.1511 


-0.4837 


10.2 


9.2 


5.8 


2 


4 


1.3 


2.3 


II 


E17 


15.1971 


-0.6266 


7.0 


5.7 




5 


10 


2.6 


3.5 


Amb. 


E18 


15.2039 


-0.6325 


8.2 


7.3 




3 


6 


2.1 


2.7 


II 


E19 


15.2052 


-0.6451 


9.2 






3 


9 


2.1 


3.7 


Amb. 


E20 


15.2336 


-0.6279 


9.2 


7.8 


3.2 


0.2 


5 


0.9 


2.5 


Amb. 


M17 EB 


E21 


15.3512 


-0.3525 


10.7 


10.1 


6.8 


2 


4 


1.0 


2.0 


II 


E22 


15.3555 


-0.2921 


9.8 


8.7 


2.5 


0.9 


6 


1.3 


2.4 


O/I 


E23 


15.3572 


-0.3770 


10.3 


9.4 


4.7 


0.9 


4 


1.0 


2.4 


Amb. 


E24 


15.3599 


-0.2973 


11.4 


10.5 


5.2 


0.2 


4 


0.2 


2.1 


O/I 


E25'' 


15.3649 


-0.4164 


7.2 


3.4 


< -0.4 


6 


10 


2.9 


3.4 


Amb. 


E26 


15.3708 


-0.3350 


7.1 


6.4 


2.4 


5 


6 


2.6 


3.0 


II 


E27 


15.3733 


-0.3714 


12.9 


12.5 


7.7 


0.2 


3 


-0.2 


1.6 


Amb. 


E28 


15.3886 


-0.3073 


10.9 


10.1 


5.9 


0.8 


4 


0.8 


2.0 


Amb. 


E29 


15.3990 


-0.3633 


12.9 


10.7 


7.1 


1 


5 


0.5 


2.7 


Amb. 


E30 


15.4620 


-0.7198 


10.3 


9.5 


6.2 


2 


4 


1.3 


2.2 


II 


E31 


15.4854 


-0.7175 


4.6 


3.7 


1.7 


9 


10 


3.6 


3.7 


II 


E32 


15.4878 


-0.4327 


12.8 


11.8 


6.9 


0.1 


4 


-0.1 


1.8 


Amb. 


E33 


15.4932 


-0.4162 


11.7 


10.1 


5.9 


0.4 


4 


0.6 


2.4 


II 


E34 


15.5115 


-0.9519 


9.9 


8.8 


6.7 


2 


7 


1.3 


3.3 


II 


E35 


15.5202 


-0.4190 


10.9 


10.0 


> 5.4 


0.7 


4 


0.8 


1.7 


II 


E36 


15.5221 


-0.4880 


8.2 


7.0 


2.4 


0.2 


1 


1.3 


1.6 


O/I 


E37 


15.5268 


-0.4078 


11.1 


11.3 


> 6.1 


0.9 


1 


1.0 


1.1 


O/I 


E38 


15.5281 


-0.4207 


12.5 


11.9 


> 4.6 


0.2 


4 


0.3 


1.7 


O/I 


E39 


15.5326 


-0.3903 


12.5 


11.2 


8.0 


0.4 


4 


0.1 


2.2 


II 


E40 


15.5330 


-0.4066 


10.6 


9.4 


> 5.3 


3 


6 


1.8 


2.9 


II 


E41 


15.5437 


-0.9338 


9.0 


7.3 


4.0 


3 


8 


2.1 


3.4 


Amb. 


E42 


15.5471 


-0.3938 


11.9 


11.5 


7.5 


0.3 


4 


0.1 


1.7 


II 


E43 


15.5489 


-1.0066 


8.3 


> 5.6 


< 0.2 


2 


7 


2.2 


2.7 


O/I 


E44 


15.5502 


-1.0063 


10.1 


8.0 


> 0.2 


1 


5 


1.4 


2.0 


O/I 


E45 


15.5553 


-0.4603 


7.1 


6.2 


2.5 


5 


8 


2.9 


3.4 


II 


E46 


15.5583 


-0.4646 


12.6 


11.0 


> 1.1 


0.2 


5 


0.4 


2.2 


O/I 


E47 


15.5584 


-0.4622 


6.5 


< 5.0 


1.1 


5 


10 


2.8 


3.8 


II 


E48 


15.5625 


-0.3870 


13.3 


11.0 


4.9 


0.1 


5 


0.4 


2.0 


O/I 


E49 


15.5691 


-0.5076 


12.0 


11.2 


7.4 


0.9 


3 


0.4 


1.4 


II 


E50 


15.6050 


-0.4788 


11.4 


10.4 


6.7 


0.6 


4 


0.5 


2.1 


II 


E51 


15.6092 


-0.4878 


12.0 


10.5 


5.4 


0.1 


4 


0.0 


2.3 


Amb. 


E52 


15.6385 


-0.5225 


10.1 


9.0 


6.0 


2 


3 


1.2 


1.7 


II 


E53 


15.6448 


-0.4847 


13.0 




6.1 


0.1 


4 


-0.0 


2.4 


Amb. 


E54 


15.6588 


-0.5279 


10.2 


9.4 


5.8 


3 


4 


1.7 


2.5 


II 


E55= 


15.6653 


-0.4989 


9.2 


9.5 


0.7 


7 


8 


2.6 


3.0 


O/I 


E56 


15.6918 


-0.8014 


11.2 


10.7 


8.3 


2 


3 


0.8 


1.5 


III 


E57 


15.7296 


-0.7243 


13.5 


11.4 


7.0 


0.2 


4 


0.0 


2.0 


II 


E58 


15.7408 


-0.6349 


12.0 


11.2 


8.1 


0.5 


3 


0.3 


1.7 


II 


E59 


15.7478 


-0.6914 


13.1 


10.7 


5.5 


0.1 


4 


0.2 


2.1 


Amb. 


E60 


15.7677 


-0.7310 


10.6 


9.4 


5.4 


1 


4 


1.1 


2.4 


II 


MC G15. 9-0.7 (Likely Unassoeiated) 


E61 


15.8856 


-0.6053 


9.0 


8.7 


6.0 


3 


4 


1.7 


2.2 


II 


E62 


15.9210 


-0.6134 


12.6 


11.9 


8.4 


0.3 


3 


-0.1 


1.6 


II 



Note. — The sources presented in this table are found within the white boundary lines but outside the 
rectangle in Fig. IIOI While the association of any given source with M17 is uncertain, the groupings of 
candidate YSOs are significant. The "Bridge" is a molecular gas structure connecting the M17 molecular 
cloud with the GMC complex extending 85 pc to the southwest. Although our cluster-finding algorithm 
selects 2 candidate YSOs apparently located in MC 15.9-0.7, these sources are likely unassociatcd. 
^ "Amb." means that the well-fit models are divided between different evolutionary Stages 
^ Source E25 is resolved by IRAC and is also an MSX point source. Its SED exhibits PAH emission. 
^ Source E55 is the candidate protostar shown in Fig. |12| 
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TABLE 3 

Other Candidate Clustered YSOs in the M17 Target Field 



Index 


I 


b 


[4.5] 


[8.0] 


[24] 




No. 


(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


Groups 


Fl 


14.5865 


-0.5708 


10.6 


9.1 


4.9 


SW 


F2 


14.5898 


-0.5822 


12.8 


11.6 


7.3 


sw 


F3 


14.5999 


-0.5422 


10.2 


9.7 


6.3 


SW 


F4 


14.6005 


-0.6248 


9.9 


9.9 


5.6 


sw 


F5 


14.6017 


-0.6177 


10.9 


10.1 


5.3 


sw 


F6 


14.6043 


-0.6229 


10.3 


9.7 


6.4 


sw 


F7 


14.6116 


-0.5520 


11.2 


9.9 


5.7 


sw 


F8 


14.6117 


-0.6259 


9.7 


9.3 


6.4 


sw 


F9 


14.6148 


-0.5843 


11.4 


10.4 


> 5.7 


sw 


FIO 


14.6235 


-0.5775 


10.7 


9.9 


> 5.5 


sw 


Fll 


14.6285 


-0.5781 


9.5 


8.6 


> 4.6 


sw 


F12 


14.6290 


-0.7614 


12.5 




5.3 


sw 


F13 


14.6297 


-0.5724 


11.1 


10.1 


> 5.4 


sw 


F14 


14.6311 


-0.5784 


10.6 


9.3 


> 2.1 


sw 


F15 


14.6314 


-0.5774 


11.7 


10.2 


2.1 


sw 


F16 


14.6327 


-0.5719 


10.9 


10.2 


> 6.4 


sw 


F17 


14.6348 


-0.8730 


13.1 


12.1 


8.1 


sw 


F18 


14.6395 


-0.7286 


11.6 


11.3 


8.3 


sw 


F19 


14.6416 


-0.8928 


13.1 


11.8 


6.1 


sw 


F20 


14.6426 


-0.7639 


11.2 


9.7 


5.3 


sw 


F21 


14.6428 


-0.8833 


12.3 


11.5 


7.2 


sw 


F22 


14.6437 


-0.9011 


10.9 


10.3 


6.8 


sw 


F23 


14.6533 


-0.7411 


12.6 


11.5 


8.2 


sw 


F24 


14.6660 


-0.8838 


12.6 


11.9 


8.0 


sw 


F25 


14.6745 


-0.8863 


10.9 


10.6 


7.0 


sw 


F26 


14.6771 


-0.7672 


9.7 


9.0 


5.1 


sw 


F27 


14.6802 


-0.8755 


11.8 


10.6 


6.0 


sw 


F28 


14.6990 


-0.8916 


12.3 


11.9 


6.9 


sw 


F29 


14.7720 


-0.3570 




8.6 


3.1 




F30 


14.7868 


-0.2176 


11.9 


10.2 


> 6.4 




F31 


14.7878 


-0.2353 


10.9 


10.5 


7.3 




F32 


14.7952 


-0.2062 


12.0 


11.7 


7.7 




F33 


14.7999 


-0.2072 


12.6 


11.1 


7.7 




F34 


14.8122 


-0.2346 


9.9 


9.3 


6.4 




F35 


14.8209 


-1.0206 


8.0 


7.4 


4.7 


sw 


F36 


14.8247 


-0.2663 


11.7 


10.2 


5.9 




F37 


14.8302 


-0.1774 


11.1 


10.7 


2.2 




F38 


14.8443 


-1.0049 


6.9 


6.4 


3.3 


sw 


F39 


14.8516 


-0.9890 


> 10.7 


10.1 


3.2 


sw 


F40 


14.8519 


-0.9927 


> 10.8 


10.0 


3.2 


sw 


F41 


14.8521 


-1.0007 


11.4 


10.7 


7.7 


sw 


F42 


14.8522 


-0.3314 


11.5 


11.0 


7.8 




F43 


14.8566 


-0.1674 


13.0 


12.1 


> 7.0 




F44 


14.8607 


-0.4358 


10.1 


9.4 


6.8 


MC G14. 


F45 


14.8621 


-0.1594 


13.8 


12.1 


7.6 




F46 


14.8624 


-0.9974 


12.8 


11.8 


8.0 


sw 


F47 


14.8637 


-0.1750 


11.0 


10.7 


7.0 




F48 


14.8695 


-0.3952 


11.3 


9.9 


5.9 


MC G14. 


F49 


14.8709 


-0.3886 


11.9 




5.5 


MC G14. 


F50 


14.8747 


-0.4154 


10.8 


9.2 


5.7 


MC G14. 


F51 


14.8840 


-0.3558 


12.4 




7.0 




F52 


14.8844 


-0.3114 


11.7 




6.3 




F53'' 


14.8891 


-0.4018 


6.1 


2.5 


< -1.0 


MC G14. 


F54 


14.8892 


-0.1977 


11.8 


11.1 


7.4 




F55 


14.8951 


-0.4004 


11.2 


7.6 


> 2.9 


MC G14. 


F56 


14.8993 


-0.2249 


10.0 


9.7 


5.2 




F57 


14.9224 


-0.2113 


12.9 


12.0 


8.1 




F58 


14.9234 


-0.2147 


12.3 




5.4 




F59 


14.9251 


-0.4058 


9.2 


8.9 


5.8 


MC G14. 


F60 


14.9261 


-0.1578 


12.4 




7.9 




F61 


14.9314 


-0.2239 


12.4 




8.0 




F62 


14.9340 


-0.1358 


11.1 


10.3 


7.0 




F63 


14.9346 


-0.2295 


12.0 


11.7 


7.7 




F64 


14.9372 


-0.4296 


10.5 


9.8 


6.2 


MC G14. 


F65 


14.9378 


-0.4229 


11.8 




7.9 


MC G14. 


F66 


14.9402 


-0.1474 


11.7 


11.1 


7.8 




F67 


14.9607 


-0.1868 


6.8 


6.1 


3.3 




F68 


14.9622 


-0.1400 


9.7 


9.4 


6.6 




F69 


14.9685 


-0.1384 


11.6 


10.8 


7.8 




F70 


14.9764 


-0.2147 


12.0 


10.5 


6.3 




F71 


14.9804 


-0.4319 


10.0 


9.2 


6.2 




F72 


14.9960 


-0.3090 


11.2 


10.8 


7.8 




F73 


15.0626 


-0.2450 


10.7 


10.2 


7.2 





TABLE 3 

— Continued 



Index 




I 


b 


[4.5] 


[8.0] 


[24] 








No. 


(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


Group'' 


F74 


15, 


.0634 


-0.3109 


10.8 


10.2 


7.4 








F75 


15, 


.0768 


-0.2528 


11.4 


10.7 


6.7 








r to 


15, 


.0773 


-U.zoii 


11. i 


1 n o 


7 1 
1.1 








Tr'zv 
r 1 t 


15, 


.0784 


-U.zo to 


1 O /I 


1 1 

11. o 


7 1 
(.1 








F78 


15, 


.0795 


-0.2609 


9.2 


8.3 


4.2 








F79 


15, 


.0817 


-0.2699 


13.0 


11.7 


> 5.2 








r ou 


15, 


.0832 


n 9fiSf^ 
-u.zuou 




Q 9 


^ 9 








F81 


15, 


.0841 


-0.1548 


10 8 


10.0 


6.5 








F82 


15, 


.1135 


n 1 

-U. 




7 R 










-T OO 


15, 


.1512 


-0 1216 


1 1 Q 


11 2 


ft ^ 








F84 


15, 


.1530 


n 1 QOQ 


12 4 


Q 7 


4 Q 








r ou 


15, 


.1684 


-U. -LZU t 


11 2 


11 4 


4 1 








r 86 


15, 


.1761 


-0.1585 


12.0 


9.4 


3.7 








F87'' 


15, 


.1831 


-0.1625 


7.8 


4.7 


< -0.5 








F88 


15, 


.1841 


-0.1589 


8.9 


6.6 


> 1.2 








F89'' 


15, 


.1850 


-0.1561 


9.7 


6.1 


> 2.3 








F90 


15, 


.1898 


-0.1400 


12.2 


11.6 


7.1 








F91 


15, 


.1904 


-0.1692 


12.6 


9.9 


> 5.9 








F92 


15, 


.2199 


-0.1451 


10.8 


9.7 


> 4.1 








F93 


15, 


.2354 


-0.1750 


9.1 


8.3 


5.7 








F94 


15, 


.5750 


-0.3101 


11.0 


10.4 


> 6.5 


Bu 


G15, 


,7 


F95 


15, 


.5754 


-0.3218 


12.3 


11.1 


7.2 


Bu 


G15, 


,7 


F96 


15, 


.5815 


-0.2963 


11.1 


9.9 


> 6.3 


Bu 


G15, 


,7 


F97 


15, 


.5869 


-0.3003 


9.1 


8.3 


3.7 


Bu 


G15, 


,7 


F98 


15, 


.5920 


-0.2923 


13.4 


10.9 


6.6 


Bu 


G15, 


,7 


F99 


15, 


.5998 


-0.2861 


10.8 


10.3 


6.6 


Bu 


G15, 


,7 


FlOO 


15, 


.6028 


-0.2040 


9.5 


8.6 


5.0 


Bu 


G15, 


,7 


FlOl 


15, 


.6051 


-0.2537 


11.3 


10.6 


6.6 


Bu 


G15, 


,7 


F102 


15, 


.6337 


-0.2282 


10.3 


9.7 


> 6.0 


Bu 


G15, 


,7 


F103 


15, 


.6408 


-0.2189 


7.9 


> 6.2 


1.0 


Bu 


G15, 


,7 


F104 


15, 


.6541 


-0.2249 


11.6 


7.9 


2.1 


Bu 


G15, 


,7 


F105 


15, 


.6603 


-0.1189 


11.4 


9.1 


5.6 








F106 


15, 


.6638 


-0.1231 


9.9 


9.1 


5.7 








F107 


15, 


.6747 


-0.1476 


10.5 




5.4 








F108 


15, 


.6768 


-0.1465 


11.9 


8.2 


4.6 








FlOO 


15, 


.6781 


-0.2651 


10.3 


9.9 


6.2 


Bu 


G15, 


,7 


FllO 


15, 


.6807 


-0.2613 


10.7 


8.9 


5.2 


Bu 


G15, 


,7 


Fill 


15, 


.6829 


-0.1424 


10.3 


9.3 


5.6 








F112 


15, 


.6929 


-0.2576 


11.4 


11.1 


7.6 


Bu 


G15, 


,7 


F113 


15, 


.6929 


-0.1355 


10.5 


7.0 


2.2 








F114 


15, 


.6956 


-0.1367 


10.8 


10.5 


> 2.2 








F115 


15, 


.7181 


-0.1648 


10.6 


9.3 


> 3.0 








F116 


15, 


.7237 


-0.2617 


11.5 


10.5 


7.3 


Bu 


G15, 


,7 


F117 


15, 


.7250 


-0.2480 


11.9 


11.5 


5.0 


Bu 


G15, 


,7 


F118 


15, 


.7574 


-0.2234 


12.2 


8.8 


4.7 








FllO 


15, 


.7580 


-0.1887 


8.9 


8.2 


5.8 








F120 


15, 


.7600 


-0.2195 


9.3 


8.6 


6.1 








F121 


15, 


.7826 


-0.3897 


12.1 


10.5 


6.5 








F122 


15, 


.7848 


-0.3864 


12.8 




7.4 








F123 


15, 


.7911 


-0.2355 


11.3 


10.3 


6.0 








F124 


15, 


.7951 


-0.4058 


> 11.8 


12.1 


7.0 








F125 


15, 


.7991 


-0.3927 


9.9 


8.8 


5.5 








F126 


15, 


.8012 


-0.3760 


8.8 


8.4 


5.6 








F127 


15, 


.8067 


-0.3818 


12.0 


11.2 


7.5 








F128 


15, 


.8108 


-0.2328 


11.7 


10.5 


8.0 








F129 


15, 


.8113 


-0.2394 


12.6 


11.1 


6.8 








F130 


15, 


.8367 


-0.2653 


10.9 


9.9 


6.3 








F131 


15, 


.8612 


-0.2350 


11.1 




5.7 








F132 


15, 


.9580 


-0.1866 


11.8 


10.7 


6.3 








F133 


15, 


.9750 


-0.2097 


9.5 


8.3 


3.9 








F134 


15, 


.9920 


-0.2097 


12.1 


10.6 


4.6 








F135 


15, 


.9930 


-0.1868 


10.9 


10.7 


6.4 









Note. — The sources presented in this table are found outside the 
white boundary hues in Fig. IIOI The majority of these sources are at 
unknown distances. 

^ Groups of these candidate YSOs are apparently associated with promi- 
nent molecular cloud structures: "SW" refers to the near end of the large 
GMC complex at — 20 km extending southwest from M17; "MC 
G14.9" refers to the molecular cloud G14.9-0.4 at — 60 km s'"^ (see 
Fig-0; and "Bu G15.7" refers to an IR/CO shell structure, 4' in diam- 
eter, centered at (/, b) ^ (15.67°, -0.29°) and u 45 km s~^. 
^ These 4 sources are resolved by IRAC. F53, F87, and F89 are detected 
by MSX and may be compact H II regions ionized by young clusters 
rather than single YSOs. 
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larger distance this source is likely to be a compact H II 
region. This interpretation is supported by the fact that 
the source also appears to have associate d radio emissio n 
from the 11-cm continuum survey of Reic h et al.l (119901 ). 
For consistency, we leave this source in the final sample 
of candidate YSOs, but we exclude it, along with other 
sources listed in Table |3] that are found outside the M17 
subregion defined in §4.1, from the analysis of YSOs as- 
sociated with M17. 

3.3. YSO Evolutionary Stages 
RW06 divided the YSO models into evolutionary 
Stages according to the envelope accretion rate Mcnv and 
the circumstellar disk mass Mdisk, both normalized by 
the mass of the central star M*: 

Stage 0/1: Mo„v/M* > IQ-^ yf-i 

Stage II: Mc„v/M, < lO'^ yr'^; A/disk/M, > 10-^ 

Stage III: M^n./M* < 10-^ yr"!; AU^^/AU < IQ-^. 

These YSO classifications are based upon model param- 
eters corresponding to physical properties, but they par- 
allel the popular observational T Tau ri classifica tion sys- 
tem that uses IR spectral indices (Ladal [T"987l) . Stage 
O/I YSOs are protostars heavily embedded in the in- 
falling gas and dust of their natal envelopes. Stage II 
YSOs, like Class II T Tauri stars, have optically thick 
circumstellar disks that dominate their near- to mid-lR 
SEDs. Stage 111 YSOs have optically-thin remnant disks 
and their SEDs are dominated by photospheric emission; 
these objects are difficult to identify by IR excess alone 
(given our error bars of > 10%-15%), hence very few ex- 
amples are found in our conservative YSO sample. Be- 
cause each candidate YSO was fit with multiple YSO 
models, we determined the evolutionary Stage of each 
source statistically. We defined the set of "well-fit" mod- 
els to each source as those models with fits meeting the 
criterion of 



-^data -^data 

where Xmin ^^i^ goodness-of-fit parameter for the best- 
fit model. We estimated the relative probability of each 
model in the set according to 

P(X') - e-(x^-x-")/2 (2) 

and normalized such that EP — 1. This allowed us to 
construct a probability distribution for the evolutionary 
Stage of each source from the Stages of all the well-fit 
models. The most probable Stage of each source was 
defined by EP(Stage)> 0.67; if this condition could not 
be met, then the Stage of the source was considered to 
be ambiguous. 

The 406 candidate YSOs in the final sample of the Ml 7 
target field are overplotted as circles on a GLIMPSE [8.0] 
image in Figure [3l The color-coding of the circles corre- 
sponds to the most probable evolutionary Stage of each 
source. Candidate YSOs are distributed throughout the 
image (the absence of sources from the top and right mar- 
gins of the image is due simply to the boundaries of our 
search area). Numerous small "clusters" of 5-10 YSOs 
are apparent, and larger groups are associated with the 
M17 molecular cloud and the rim of M17 EB. Candidate 



YSOs associated with the M17 molecular cloud would 
dominate the sample if not for the strong masking effect 
of the bright diffuse mid-lR emission from the H II region 
and PDR. The point-source detection limit rises dramat- 
ically near the M17 H II region, and only the most lumi- 
nous YSOs are detectable within the area where the 8.0 
image in Figure [3] appears black. The MIPSGAL 24 
fim image is completely saturated within the H II region, 
and a relatively high fraction of candidate YSOs near 
the H II region have ambiguous Stage (magenta circles) 
because 24 /im fluxes are an important discriminator be- 
tween disk- and envelope-dominated sources (RW06). 

3.4. Mid-IR Color-Color Plots 

In Figures [4] and [5] we compare our model-based clas- 
sifications of all stars and YSOs in the M17 target field 
with classification systems derived from IRAC mid-IR 
colors (Allen et al. 2004, hereafter A04; Gutermuth et 
al. 2008, hereafter G08). A04 showed, through a study 
of several nearby young clusters, that a box-shaped re- 
gion in the IRAC [3.6] — [4.5] versus [5.8] — [8.0] color-color 
diagram (shown in Figs. [4|) contains primarily Class II 
YSOs with disks. A04 found that protostars were lo- 
cated in the redward extension of this "disk domain" 
along both color axes. G08 built upon the work of A04 
to devise a more sophisticated system for the classifica- 
tion of YSOs in the young embedded cluster NGC 1333, 
utilizing multiple IRAC color-color planes. In Figures H] 
and [5] the YSO selection and classification criteria of A04 
and G08 are compared to our YSO selection and Stage 
classification using the RW06 models. Sources plotted in 
Figures SJa, b and ^jp., b are color-coded according to their 
most probable evolutionary Stage from the RW06 mod- 
els: red, orange, green, and lavender for Stages 0/1, 11, 
III, and ambiguous, respectively, while gray dots repre- 
sent sources best fit by stellar photospheres, and blue 
dots show sources that were well-fit by YSO models but 
were excluded from the final sample (see §3.2 above). 
Figures ^ and show only sources for which the classi- 
fication using model fitting differs from the classification 
using the G08 color-color criteria (sources with ambigu- 
ous Stage are not plotted). In Figures IDs and[51c, sources 
are plotted according to their G08 color classification: 
protostars are red while Class II sources are orange. Like 
virtually all classification systems based upon colors, nei- 
ther the A04 nor the G08 criteria allow ambiguity be- 
tween protostars and Class II sources, but we stress that 
ambiguity is an inherent property of any YSO taxonomy, 
whether based upon colors, spectral indices, or SEDs. 

Although the A04 disk domain contains many disk 
sources, it does not contain all candidate disk sources 
and selects a significant number of diskless stars with 
reddened photospheres. The locus of Stage II sources 
in Figure SJa is redder than the disk domain by [3.6] — 
[4.5] ~ 0.3 mag . The same displacement was found by 
llndebetouw et all ([2007) for YSOs in the Eagle Nebula 
(M16), another massive star formation region near the 
Galactic mid-plane, located ~2.5° from M17. The dis- 
crepancy could be largely resolved by reddening the disk 
domain by an amount equivalent to Ay ~ 20 mag, a rea- 
sonable modification considering that the regions studied 
by A04 were not highly-reddened massive star formation 
regions in the Galactic plane. The G08 color classifica- 
tion selects 1882 Class II sources and 57 protostars that 
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16.000 15.800 15.600 15.400 15.200 15.000 14.800 14.600 

Galactic longitude 

Fig. 3. — GLIMPSE 8.0 image of the M17 analysis target field (inverted logarithmic grayscale). Positions of the 406 
candidate YSOs in the final sample are marked by circles color-coded according to the most probable evolutionary Stage of 
each source: red circles are Stage O/I candidates (deeply embedded in infalling envelopes); yellow are Stage II (optically thick 
circumstellar disks); green are Stage III (optically thin circumstellar disks); and magenta are candidate YSOs for which the 
Stage is more ambiguous. Blue squares denote candidate AGB stars removed from the YSO sample. The apparent lack of 
candidate YSOs in the extreme top and right portions of the image reflects the boundaries of the search area. 



wc find to be well-fit by reddened stellar photospheres 
{gray dots in Figures H^; andO*). These sources are a 
striking feature of Figure [3]c, where they are found pre- 
dominantly to the left of the disk domain, in the locus of 
photospheres (AOS). Largely on account of these sources, 
the total number of YSOs selected by the A04 criteria or 
especially the G08 criteria is far greater than the number 
of candidate YSOs in our conservative sample. Again, 
the fact that the G08 color criteria were derived through 
the study of a nearby star formation region located away 
from the Galactic plane would explain the large num- 
ber of reddened stellar photospheres selected as YSOs. 
In addition, both the A04 and G08 criteria were derived 
from studies of regions that, unlike M17 and M16, lack 
many massive YSOs, which should contribute to the dis- 
parities between the color-colo r classifications a nd the 
model-based classifications (Whitnev et al.|[2004[ ). 

When compared to our evolutionary Stage classifica- 
tions, the GOB classification system confuses a signifi- 
cant number of protostars for more evolved YSOs, and 
vice versa. This disagreement is evident in Figure 31 but 
it can be more readily explained by examination of Fig- 
ure \S\ in which sources are plotted using a color-color 
diagram that combines IRAC with MIPS [24] photom- 
etry. As shown by RW06, this color space illustrates 
the value of 24 fim data for distinguishing between pro- 
tostars (Stage O/I) and disk-dominated (Stage II) YSOs. 
Figureinh shows a progression toward less evolved sources 
with increasing [8.0] — [24] color, and many sources with 



ambiguous Stage determinations (lavender) lie along the 
boundary between the Stage O/I and II populations. In 
Figure O?, all Stage II sources {orange dots) plotted for 
[3.6] — [5.8] > 1.5 are protostars according to the G08 
classification system (these sources are red dots in Fig. 
Oc), which relies upon photometry at wavelengths <8.0 
fim. 

Two main classes of "false" YSOs are illustrated in the 
IRAC color-color plots of FiguresHl (1) Candidate AGB 
stars (open blue circles) are intermingled with YSOs and 
appear very difficult if not impossible to distinguish by 
IRAC colors alone. (2) The lower-right part of Figure 
Ha, centered on [3.6] - [4.5] = for 1 < [5.8] - [8.0] < 2.5, 
contains numerous sources well- fit by YSO models {blue 
dots) that were discarded from the final sample because 
the IR excess emission appears at [8.0] only (see §3.1 
above) . 

Spurious correlations of 24 /im sources with IRAC 
sources can also create artificial IR excesses (see §3.1), 
and consequently many sources discarded from our fi- 
nal YSO sample appear in Figure [SJa as blue dots near 
[3.6] - [5.8] = for 2 < [8.0] - [24] < 4. Candidate AGB 
stars are better distinguished from YSOs in Figures [S] 
than in Figures [H Three candidate AGB stars (see §3.2 
above) with [3.6] - [5.8] > 1.5 but [8.0] - [24] < 1 are 
located far from YSOs and stellar photospheres. Sev- 
eral other candidate AGB stars have [8.0] — [24] ~ 2, 
placing them on the blue edge of the Stage II YSO lo- 
cus. This supports the interpretation that AGB stars 
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Fig. 4. — Mid-IR color-color space showing sources in the 
target field that were detected in all 4 IRAC bands. The "disk 
domain" of A04 is contained within the black box. (a) Can- 
didate YSO selection and classification from model fitting. 
Small gray dots are sources well-fit by stellar photosphere 
SEDs. Large dots are sources in the final YSO sample {red: 
Stage O/I; orange: Stage II; green: Stage III; and lavender: 
ambiguous), (b) Sources for which our model-based classifi- 
cation differs from the IRAC color-color classification of G08. 
Sources discarded from the final YSO sample (small blue dots) 
remain in this plot if they have IRAC colors consistent with 
YSOs according to G08. (c) G08 YSO selection and classifica- 
tion (gray: stellar photospheres; red: protostars; and orange: 
Class II). In all plots, blue circles are candidate AGB stars, 
reddeni ng vectors for Ay = 30 mag based on the ext inction 
laws of IWeingartner fc Draind (|2001l ) and Undcbctou w et al.l 
are shown as open and filled arrows, respectively, and 
black crosses give typical photometric errors. 
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Fig. 5. — Same as Fig. |4j except a color space combining 
IRAC with MIPS [24] photometry is shown. 



are typically surrounded by dust shells with smaller ef- 
fective radii than the disks and envelopes of YSOs (see 
iRobitaille et al.ll2008l . for example). 

We have experimented with the A04 and G08 color- 
color criteria for selecting and classifying candidate YSOs 
because these are among the most well-developed anal- 
ysis techniques available for IRAC. It is not appropri- 
ate to apply these precise color cuts directly to the in- 
terpretation of sources in the M17 region, because they 
were developed through the study of relatively nearby 
{d < 1 kpc) regions that, unlike M17, are located in 
relatively low-density molecular clouds away from the 
Galactic plane. This illustrates one of the principal ad- 
vantages of the model- fitting approach. All methods for 
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Fig. 6.— GLIMPSE 8.0 fim image of the control field, with 
the same color overlays as in Fig. [3] The area searched for 
YSOs is shown by the black circle. The control subregion 
chosen to mirror the position of M17 about the Galactic mid- 
plane is the boxed area, 0.5° on a side. 

selecting populations of objects based on broadband photo- 
metric colors are based upon models. The RW07 model 
fitter, with its integrated capacity to apply interstellar 
extinction to SED models before matching them with ob- 
served fluxes, can be employed consistently for the study 
of star formation in various environments without the 
need to empirically redraw the selection criteria. A sec- 
ond major advantage is the ability to use all available 
photometric information simultaneously to constrain the 
physical parameters of large numbers of sources. For 
mo re on the advantages of the model-fitting approach, 
see llndebetouw et all ()2007t ). 

3.5. Contaminants in the YSO Sample 

Even the most carefully selected sample of candidate 
YSOs will identify other objects that have similar IR 
colors. External, dusty galaxies are the most common 
contaminant in YSO searches using deep IRAC obser- 
vations of nearby star formation regions a way from the 
Gala ctic plane, such as the c2d Survey (jPorras et al.l 
[2nnl. GLIMPS E is a shallow survey, and employing 
the IPorras et al.l ((2007) criterion for galaxy selection to 
our final M17 target sample of YSOs yields no galaxy 
candidates (the control sample contains a single source 
that may be a galaxy). In the Galactic plane, the objects 
most likely to masquerade as YSOs in significant num- 
bers are evolved stars, specifically AGB stars with dusty 
winds and unresolved planetary nebulae (PNe). In ad- 
dition, because we are primarily interested in analyzing 
the population of YSOs associated with M17, YSOs in 
the foreground or background must also be treated as 
contaminants to our sample. 

A GLIMPSE 8.0 ^m image of the control field is pre- 
sented in FigurelU with the final sample of 106 candidate 
YSOs and AGB stars overplotted as in Figured Can- 
didate YSOs are distributed throughout the search area 
with an average density of 140 sources deg~^, or nearly 
60% of the average source density in the target sample 
(see Table [IJ. The source density in the control sample is 
not uniform, decreasing with increasing Galactic latitude 



across the 1-degree span of the field, which implies a low 
Galactic scaleheight for the source population. Because 
of this variation, we define a box 0.5° on a side centered 
at {l,b) — (15.1,0.7) in the control field, to sample di- 
rectly the density of contaminants at Galactic latitudes 
mirroring the location of M17. The source density in this 
control subfield is 100 deg-^ (Table [T]). 

The control sample represents the contamination ex- 
pected to be present in the sample of candidate YSOs 
in the M17 target field. To estimate the numbers of 
evolved stars expected in the control and target samples, 
we turned to a detailed model of the IR point source sky 
(|Wainscoat et al.l [1991 : ICohenI [l993l) . This Sky model 
incorporates 87 different source populations, including 
PNe and several types of AGB stars, into a geometric 
representation of the Galactic bulge, halo, and thin and 
thick disks. The model has been shown to a gree well 
with /i?^S' observations (jWainscoat et al.|[l992[ ). 

Cumulative source densities as functions of [4.5] mag- 
nitude are presented in Figure [71 Curves generated by 
the Sky model for oxygen- rich AGB stars (AGB M) , op- 
tically visible/obscured carbon stars (AGB C/CI), and 
PNe are plotted along with YSO candidates detected in 
the target and control samples and subsamples (see §4). 
One striking feature of these plots is the high bright- 
ness of AGB stars in the IR. The brightest magnitude 
plotted in Figure [7| is magnitude 8, because nonlinear- 
ity begins to affect the [4.5] band for brighter sources. 
The source density of all AGB types brighter than mag- 
nitude 8 is more than twice the density of comparably 
bright sources selected for our YSO sample, even before 
suspected AGB candidates are removed (dashed curves 
in Fig. [7]). While the brightest AGB stars will saturate 
the IRAC detectors and hence not be extracted into the 
point-source Archive, there are too few such sources vis- 
ible in the target and control fields to account for a sig- 
nificant number of the AGB stars predicted by the Sky 
model. While most AGB stars in the inner Galaxy will 
be in the GLIMPSE Archive, evidently not all types of 
AGB star have large IR excesses that allow them to be 
confused with YSOs. 

Removing suspected AGB stars from the control field 
sample cut the cumulative source density at magnitude 
8 by more than half (the difference between the dashed 
and solid red curves in the right panel of Figure [7]) . As 
the cumulative source density with magnitude plotted in 
Figure [7| shows, the Sky model predicts that <30% = 
(70 deg-2 - 50 deg"^)/70 deg"^ of 0-rich AGB stars 
{blue curves) and virtually no carbon stars {green curves) 
are fainter than [4.5] = 8 mag. Suppose that all the IR 
excess sources remaining in the control sample brighter 
than [4.5] = 8 mag are 0-rich AGB stars, with no fainter 
AGB stars having been removed. Scaling down the black 
curve to match the solid red curve at 8 mag in the control 
field plot of Figure [7| shows that the final control sample 
could still contain up to 15 AGB stars deg^^, or ^12% 
of the overall source density (~125 deg~^; value of the 
solid red curve at [4.5] — 14 mag). This represents the 
maximum possible fraction of AGB stars remaining in 
the control sample of candidate YSOs. The fraction of 
AGB stars remaining in the M17 target sample will be 
even lower (~7%), because the overall source density is 
1.7 times higher (Table [T]). 
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Fig. 7. — Cumulative densities of different types of red sources as a function of [4.5] magnitude. The curves for AGB stars 
(AGB M are 0-rich and AGB C fc CI are optically visi ble & obscured C-rich AGB stars, respectively) and PNe were generated 
by the Sky model (jWainscoat et al.lll992l : ICohenlllQOaD . The sohd YSO curves represent the target and control samples and 
subsamples, while the dashed curves show the solid curves before bright candidate AGB stars were removed. 



There remains some potential contamination from un- 
resolved PNe. According to the predictions of the Sky 
model (Fig. [7]), these represent < 10% of sources in the 
control sample and < 5% of sources in the target sample. 

We conclude that the majority (>80%) of sources se- 
lected as candidate YSOs in the control field are ac- 
tual YSOs, albeit at unknown distances. This is higher 
than the 60% fraction of YSOs found in the recent 
cat alog of dusty r e d sou rces in GLIMPSE compiled 
bv iRobitaille et aP ()2008l ). This catalog contains only 
sources with [8.0] > 10 mJy. Virtually all Galactic AGB 
stars meet this brightness criterion, but only ~50% of the 
sources in our control sample do. Th erefore, if we apply 
the [8.0] cut of lRobitaille et al.l ()2008f ) to our control sam- 
ple, the fraction of candidate evolved stars doubles, from 
20% to 40%, in good agreement with the new red source 
catalog despite the small size of our control sample. 

The density of IR excess sources drops noticeably 
across the 1° increase in Galactic latitude on either side 
of the mid-plane (Figs. [3] and [6l). This suggests that star 
formation, in the form of isolated YSOs or small clus- 
ters, is distributed throughout the Galactic plane with a 
low scaleheight. The number of candidate YSOs in the 
M17 vicinity represents a significant enhancement over 
this Galactic background YSO population, as we discuss 
in the following section. 

4. THE EXTENDED YSO POPULATION AND 
MOLECULAR CLOUD STRUCTURE OF M17 

A ^^CO zero moment map for v = 12-26 km s~^ and 
a ^^CO channel map at w = 19 km s~^ are compared to 
the 8 /Ltm image of the M17 target field in Figure [H The 
^^CO and ^^CO maps from the HHT reveal the underly- 
ing structure of the molecular gas outlining the extended 
regions of M17 along with M17 EB and portions of MC 
G15.9-0.7. M17 EB is a coherent, albeit clumpy, kine- 
matic structure in molecular gas at w = 19 km s^^. The 
CO data thus confirm that the entire dust shell observed 



in the IR has a common velocity with M17, supporting 
our conclusion that M17, M17 EB, and MC G15.9-0.7 
are all part of the same dynamical complex. M17 EB, as 
a bubble, should exhibit a 3-dimensional, roughly spher- 
ical structure. Many astrophysical shell structures ex- 
hibit a ring-like morphology in optically thin emission 
due to limb-brightening, and emission from the interior 
of the ring may be difficult to detect. Integrating the 
^•^CO data cube spatially over a 17' x 11' box centered at 
{I, b) = (15.6°, —0.7°), we definitively detect an emission 
line near v = 20 km s~^ with a peak at (Is) = 0.05 K. 
The line does not exhibit the double-peaked profile of an 
expansion signature, which means that if the bubble is 
still expanding, the expansion velocity is low, less than a 
few km s~^, and the line profile is dominated by internal 
motions within the shell. 

4.1. Correction for Contamination from Unassociated 

Sources 

We have claimed that the control field sample of can- 
didate YSOs mirrors the population of unassociated 
sources contaminating the M17 target field sample. We 
can use the control sample to remove contaminants sta- 
tistically from the target sample, but first we must ver- 
ify this claim. The molecular cloud complex associated 
with M17 contributes a significant enhancement to the 
gas column density in the target field with respect to 
the control field. This additional source of interstellar 
extinction could bias the source selection between our 
flux-limited samples. 

4.1.1. Extinction Map 

To estimate the extinction through the molecular 
cloud, we employ a model for CO line excitation based 
on a n escape probability formalism for the CO line radia- 
tion (|Kulesa et al.l [2005). The observed ^^CO and ^^CO 
( J = 2 — > 1) line intensities, velocity widths, and peak 
brightness temperatures constrain the model at each 
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Fig. 8. — CO maps of extended molecular cloud structures associated with M17. {Top) HHT ^^CO moment map for v = 12-26 
km s~^ The image colorbar is in units of K(r4) km s ^. Contours show 90 cm continuum observations from the VLA (at 
intervals of 5, 10, 20, 30, 40, 60, and 80% of the map peak) that trace thermal emission from ionized gas associated with the 
M17 H II region. {Bottom) GLIMPSE 8 ^^ra image. Contours show HHT ^^CO emission at u = 19 km s"^ 
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map pixel to give the total CO column density, A^(CO), 
and the mean gas kinetic temperature, {Tkin)- We es- 
timate the mean H2 density, (n(H2)) from the derived 
A^(^'^CO) column density. (n(H2)) varies from 10^ cm~'^ 
in the extended cloud structures to 10^ cm~^ in the M17 
South core near the H II region. The derived A^(CO), 
(n(Il2)), and (Tkin) are th en applied as constraiii t s on a 
CO PDR model based on Ivan Dishoeck fc Blackl ()1988D 
to estimate the far UV radiation field, Go, incident on 
the molecular gas, and the total hydrogen column den- 
sity 7V(H) = A^(H'') + 2A^(H2) (H-atoms cm^^). The 
corresponding extinction Ay, is calculated from iV(H) 
using the value N(}i)/Ay = 1.9 x 10^^ cm~^ mag~^ 
(iBohhn. Savage, fc Drakj|1978[ ). Our analysis thus at- 
tempts to incorporate the variations in both gas density 
and the far UV radiation throughout the cloud complex. 
The resulting distribution of extinction is shown in Fig- 
ure 9, as a map of Ay over the region mapped in ^^CO 
and "CO. 

The PDR model outputs the total hydrogen column 
density used in the extinction calculation. Abundances 
measured from the M17 H II region give X — 0.698 
(jCarigi fc Peimber3l2008l ). and including unobserved He 
and metals the entire region mapped by HHT contains 
1.35 X 10^ Mq total gas mass in the velocity range occu- 
pied by M17. 

The extinction map in FigurelHlrepresents the differen- 
tial extinction between the target and control fields due 
to the molecular and atomic gas associated with M17. 
The molecular gas occupies less than half of the image 
area, while candidate YSOs are found throughout (Fig. 
[3]). Much of the M17 molecular cloud has Ay > 15 mag, 
reaching Ay > 90 mag in the dense core of M17 South. 
The extended outer regions of M17 and MC G15.9-0.7, 
however, produce Ay < 15 mag. This correspo nds to 
Ak < 2 mag and ^[4 5] < 1 mag (Indcbctou w et al.l 
|2005() generally, and much less through the interior of 
M17 EB and other areas lacking high-density molecular 
gas. The map in Figure [H represents the spatial average 
of Ay over each beam. Hence the extinction could be 
higher along localized sightlines through compact, dense 
clumps of molecular gas that are smaller than the 32" 
(FWHM) resolution of our map (0.33 pc at 2.1 kpc), 
but such clumps occupy an insignificant fraction of the 
mapped area. We thus do not expect the differential ex- 
tinction between the target and control fields to produce 
a significant disparity in mid-IR point source detections, 
except in the M17 molecular cloud itself. 

The Ml 7 molecular cloud and PDR do create signifi- 
cant, spatially varying extinction and, more importantly, 
produce bright diffuse mid-IR emission that saturated 
the MIPSGAL 24 /im images and drastically reduced the 
point-source sensitivity of GLIMPSE. The surface den- 
sity of GLIMPSE Archive sources is 14% higher in the 
control field than in the target field. Much of this dis- 
crepancy is due to the sharp drop in source detections 
near the M17 H II region and molecular cloud, which oc- 
cupy '-^10% of the area in the target field. We therefore 
consider the control field ill-suited for comparison with 
the central regions of M17, but for the remaining 90% 
of our extended target field the control sample gives an 
accurate picture of the contamination to the target YSO 
sample. 



4.1.2. Selection of YSOs Associated With M17 

The prevalence of candidate YSOs distributed 
throughout the M17 target field requires that a distinc- 
tion be drawn between areas that are likely dominated by 
sources associated with M17 and external star formation 
regions. Using the distribution of 8 /im emission along 
with the CO velocity structures revealed by our HHT 
data, we defined a subregion of the target field that is 
dominated by molecular gas associated with M17 (area 
bounded by white in Fig. [T0| . The 195 candidate YSOs 
within this region make up the M17 target subsample 
(Table [ij. This subsample is directly comparable to the 
control subsample found within the boxed region in Fig- 
ure [6] 

In Figure [TOl the background contamination to the 
YSO sample has been statistically subtracted such that 
only candidate YSOs exhibiting a significant degree of 
clustering are plotted. This was accomplished using 
a simple "friends-of-friends" nearest-neighbors analysis. 
We first computed the median angular separation ^con 
between each source in the control field sample and its 
N nearest neighbors, where 3 < TV < 10. We then cal- 
culated ^targ, the corresponding nearest-neighbors sep- 
aration for each of the candidate YSOs in the target 
sample. The median value of ^con for all sources in the 
control sample is Ocon- Finally, we selected the subset 
of sources in the target field for which ^targ < /©con- 
The factor / comes from the ratio of surface densities 
between the target subsample and the control subsample 
(/ = [Scon/Starg]°'^), with / = 0.7±0.1 usiug the values 
from the bottom row of Table [1] and assuming counting 
statistics. We found that the subset of candidate YSOs 
selected by this procedure was generally insensitive to the 
choice of N for > 3, while the uncertainty in / cre- 
ates a ~20% uncertainty on the total number of sources 
selected. The final sample of candidate YSOs shown out- 
side the white rectangle in Figure [TO] was selected using 
N = 5, for which Ocon = 4.28'. For the central parts 
of M17 we did not correct for contaminants because our 
control field is unsuitable, as noted above. The 34 can- 
didate YSOs found within the white rectangle in Figure 
[To] are presented in Tabic [4] Excluding these 34 sources, 
we have discarded 197/372 = 0.47 of the sources from 
the original target sample. Given that we expect a frac- 
tion p ^ 100/205 = 0.49 of the sources in the target 
sample to be contaminants (Table [1]) , our cluster-finding 
procedure has succeeded in removing the expected con- 
tamination. 

After correcting for contaminants the final M17 YSO 
subsample contains 96 sources (Tables [2] fcHJ). Two main 
concentrations of candidate YSOs within the final sub- 
sample are apparently associated with M17 EB and the 
M17 H II region, respectively (Fig. [TO]) . The final sub- 
sample excludes 135 candidate YSOs (Table [3]) that ei- 
ther are unassociated with M17 or are associated with 
the extended molecular c l oud complex to the southwest 
(lElmegreen fc Ladal[l976l : iSanders et al.l[l98^ . YSOs in 
the latter group are at the same distance as M17 and are 
associated with the M17 complex, but they belong to an 
additional extended YSO population that is beyond the 
scope of this work. 

4.2. Star Formation in Mil EB and MC G15. 9-0.7 
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4.2.1. YSO Mass Function 

We constructed a mass function for the observed YSO 
population following a proc edure employed previo usly by 
IShepherd et~al] (|2007D and lWhitnev et all pOOSD . Using 
Equation 2, we assigned a x^-weighted probability for 
the stellar mass of each well-fit model to a given 
candidate YSO. By summing over the probabilities of 
all well- fit models (defined by Equation 1) and normal- 
izing to unity, we derived a probability distribution of 
Mi, for each candidate YSO. This procedure is similar 
to that employed to find the probability distributions of 
evolutionary Stage (§3.1), but a finer binning was used. 
Results of fitting models to 62 sources in the extended 
Ml 7 region (within the white boundary lines in Figure 
[TUl but excluding the white rectangle) are presented in 
Table [21 Ranges giving 95% confidence intervals for M^, 
and Ltot (total bolometric luminosity, including lumi- 
nosity produced by accretion) for each candidate YSO 
are listed in Table [21 Summing the probability distribu- 
tions of over the 62 sources in Table [H produces the 
YSO mass function (YMF) plotted in Figure [TTI 

The YMF in Figure ITT] exhi bits . for intermediate-mass 
sources, the ubiquitous iSalpeteii (1955) power-law slope. 
The turnover for Mj, < 3 Mq is not a real break in 
the YMF but instead reflects the incompleteness of our 
YSO sample with respect to the detection of lower-mass 
sources. A completeness limit of ~3 M© is reasonable 
for the detection of disk-dominated sources. The [8.0] 
flux density of an AO V photosphere w ith M^, = 3 .2 M©, 
Teff = 10,800 K, and = 2.5 R© ([Coil2000[ ) at 2.1 
kpc is 3.9 mjy, or '^1 dex above the nominal detection 
limit of the point-source Archive. The extended M17 
field from which the YMF has been sampled exhibits nei- 
ther high source confusion nor very bright diffuse 8 /xm 
emission, hence our YMF represents a nearly complete 



sample of intermediate-mass IR excess sources that will 
become main-sequence types AO or earlier. 

We estimated the size of the larger YSO population 
represented by o ur intermediate-mass YSO sample by 
scahng the IMF of lMuench et al.l (|2002f) to fit the part of 
our YMF (3 < < 8 M©) where it is most complete 
(solid curves with dash-dotted uncertainty envelopes in 
Fig. fTT]) . This procedure generated a model IMF, ^(M*), 
for the YSO population. ^{M-i,) has a power-law slope of 
F = -1.2 ± 0.2, as defined by 



Mi 



rf(logM,)' 



for M^ > 0.6 M©. After the break poi nt, C(M^) flat- 
tens t o r = -0.2 ± 0.2. Although the full iMuench et all 
(|200l IMF is a 4-part power law extending to brown 
dwarf masses, we use only ^(M^ > 0.1 M©) for consis- 
tency with the BFT07 results^'' (see §4.3 below). Inte- 
grating ^(M^, > 0.1 M©), we predict that the extended 
YSO population of M17 numbers 950 ± 90 and repre- 
sents a total stellar mass of 750 ± 50 M©. These val- 
ues are lower limits. We cannot be certain that our fi- 
nal YMF is 100% complete even for intermediate-mass 
sources. We treat any unresolved multiple sources as sin- 
gle YSOs, which means that we are generally sensitive 
only to the most luminous component of a multiple sys- 
tem, given that Ltot oc M° wi th a > 2.5 for PMS 
stars (jBernasconi Sz Maeded Il996l ). The mass deter- 
mined from the YSO fits will be slightly higher than the 
most massive source in the beam but substantiall y lower 
than the combined mass of the sources ([Whitney et all 
I2OOI . 

^■^ Extrapolating down to the hydrogen-burning limit at 0.08 
Mq increases neither the total number of YSOs nor the total stellar 
mass by a significant amount. 
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Fig. 10. — Candidate YSOs exhibiting significant clustering (see text). Circles mark YSO positions and are color-coded 
according to evolutionary Stage as in Fig. [3] The 8 /im image and CO contours are the same as in Fig. [8] (bottom). The thick 
white boundary lines enclose the CO emission and pass through breaks in the spatial distribution of YSOs. The thick white 
rectangle shows the field that has been observed by Chandra, within which bright diffuse mid-lR background emission severely 
affects the IRAC point-source sensitivity. Green rectangles show the fields enlarged in Figs. ll2lfc [T3l 
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Fig. 11. — Mass function of the 62 candidate YSOs found 
inside the white boundary lines but outside the white rectangle 
in Fig. 1101 (hi stogram with error b ars). The Orion Trapezium 
cluster IMF (jMuench et al.|[2002l 'l has been scaled to match 
the mass function for ~3 M© < M* < 8 Mq (heavy curve 
with dash-dotted curves showing uncertainty envelopes). 

The YMF in Figure [11] contains candidate YSOs dis- 
tributed throughout a large volume of space that can- 
not be regarded as members of any single cluster. The 
YMF is, however, drawn from a coeval stellar population 
comprising the most recent generation of star formation 
in the extended M17 target field (excluding YSOs as- 
sociated with the central NGC 6618 cluster). Because 



these sources were selected on the basis of mid-IR ex- 
cess emission, and a negligible fraction were fit with 
Stage III YSO models, the YMF population is charac- 
terized by optically-thick circumstellar disks and/or in- 
falling envelopes. The disk destruction timescale is a 
decreasing function of the mass of the central star, with 
solar- mass YSOs losing their disks on timescales of <2 
Myr (|Haisch. Lada. fc Ladall2001al) . The intermediate- 
mass YSOs dominating our sample are expected to evolve 
more quickly, making them even younger. 

4.2.2. Timescales of Evolutionary Stages and Lower Limit 
on Star Formation Rate 

The duration of the Class I phase for YSOs in low-mass 
star formation regions is "^O-l Myr (|Kenvon &: HartmannI 
[1991 iHaisch. Lada. k Ladal[200l . We can make an in- 
dependent estimate of the duration of the Stage O/I phase 
from the accretion ages of the well-fit models for each 
source. We define "accretion age" as Ia — 10^^A/*/Afcnv 
[Myr] . The accretion age is a proxy for the evolutionary 
age i of a YSO, assuming both that Mem is constant, on 
average, over tA and that all of the accreting gas eventu- 
ally reaches the central star. In reality, average accretion 
rates decrease with time, causing tA to overestimate t, 
but accretion-driven outflo\ ys prevent 3-30% of the in- 
falling gas from accreting ([Richer et al.l 120001 ) , driving 
tA in the opposite direction. Combining the probability 
distribution of tA for the 25 candidate Stage O/I YSOs in 
Tables[2]&[31 we find that tA has a median value of 0.075 
Myr, in good agreement with other estimates of Class I 



The Star Formation History of M17 



19 



TABLE 4 

Candidate YSOs in the Region Observed By Chandra 



Index 
No. 



I 

(dog) 



b 

(dog) 



[4.5] 
(mag) 



[8.0] 
(mag) 



(Mo) 
Min Max 



log(LTOT/L0) 
Min Max 



Evolutionary 
Stage'' 



With X-ray Counterparts 



XI 


14.9846 


-0.6718 


9.9 


8.7 


2 


3 


1.2 


1.6 


II 


X2'° 


14.9855 


0.6716 


10.6 


9.6 


1 


5 


1.0 


2.5 


II 


X3 


14.9954 


-0.6332 


8.4 


5.7 


1 


7 


2.1 


2.8 


O/I 


X4<= 


14.9962 


-0.6733 


< 3.6 




9 


10 


3.7 


3.9 


II 


X5 


14.9994 


-0.6304 


7.3 


6.9 


5 


6 


2.6 


3.2 


II 


X6 


15.0198 


-0.6768 


8.6 




4 


10 


2.3 


3.8 


II 


XT'' 


15.0386 


-0.6513 


6.2 




6 


15 


3.0 


4.4 


Amb. 


X8 


15.0470 


-0.7023 


9.9 


8.3 


1 


5 


1.1 


2.7 


Amb. 


X9 


15.0530 


-0.6934 


7.1 


6.0 


5 


9 


2.6 


3.6 


Amb. 


XIO 


15.0568 


-0.6068 


9.9 


8.1 


2 


5 


1.5 


2.0 


O/I 


Xll 


15.0699 


-0.6646 


9.4 


6.8 


1 


5 


1.3 


2.4 


O/I 


X12 


15.1146 


-0.6079 


11.4 




1 


7 


1.1 


2.6 


O/I 


X13 


15.1613 


-0.6335 


10.4 


9.5 


0.5 


4 


0.8 


1.9 


Amb. 


X14'' 


15.1772 


-0.6593 


9.0 


8.2 


2 


5 


1.5 


2.3 


II 


Without X-ray Counterpart; 


i 
















NXl'i 


14.8923 


-0.6556 


11.0 


10.1 


2 


4 


1.1 


2.0 


II 


NX2d 


14.9277 


-0.6740 


8.0 


7.5 


4 


5 


2.3 


2.7 


II 


NX3<* 


14.9278 


-0.6473 


11.5 


10.3 


0.4 


7 


0.7 


2.7 


O/I 


NX4 


14.9761 


-0.7292 


11.3 




0.4 


7 


0.8 


3.0 


O/I 


NX5<= 


14.9791 


-0, 


.6651 


6.6 


3.3 


8 


11 


3.5 


3.9 


Amb. 


NX6 


14.9806 


-0, 


.6787 


11.0 


9.9 


0.5 


5 


1.0 


2.4 


Amb. 


NX7 


15.0494 


-0, 


.6838 


7.6 




4 


10 


2.2 


3.8 


Amb. 


NX8 


15.0576 


-0, 


.6830 


10.8 


8.7 


0.6 


7 


1.0 


2.9 


O/I 


NX9 


15.0589 


-0, 


.5735 


10.8 


9.9 


0.5 


4 


0.7 


1.5 


Amb. 


NXlQf 


15.0848 


-0, 


.6631 


7.7 


6.5 


2 


10 


2.1 


3.6 


Amb. 


NXll 


15.0895 


-0, 


.5641 


11.0 


9.2 


0.7 


5 


1.0 


2.6 


O/I 


NX12 


15.0984 


-0, 


.6658 


6.3 


6.0 


5 


14 


2.9 


4.2 


III 


NX13 


15.1176 


-0, 


.5618 


9.2 


8.3 


3 


7 


2.0 


2.7 


O/I 


NX14d 


15.1337 


-0, 


.8681 


10.3 


8.5 


0.4 


6 


1.0 


3.0 


II 


NX15<* 


15.1348 


-0, 


.8683 


13.0 


10.7 


0.2 


6 


0.3 


2.3 


O/I 


NX16 


15.1366 


-0, 


.6657 




5.8 


7 


14 


3.2 


4.2 


O/I 


NX17 


15.1784 


-0, 


.6299 


10.5 


9.6 


0.3 


3 


0.7 


1.6 


Amb. 


NX18 


15.1808 


-0, 


.6222 


8.2 


6.7 


0.5 


8 


1.5 


2.9 


O/I 


NX19<* 


15.1844 


-0, 


.8984 




5.6 


5 


8 


2.9 


3.3 


II 


NX20 


15.2017 


-0, 


.8404 


8.5 


7.5 


2 


9 


2.0 


3.2 


O/I 



The "Amb.' 



same as m 



Tabic [2] 

"protostais" in Table 7 of BFT07. 



Stage designation is the 
^ Selected as X-ray emitting candidate 

X4 is the KW Object (see Fig.fTst. 
^ These 6 sources are located away from the bright 8 /im emission of the M17 PDR, increasing the odds 
that they are unassociated contaminants. 
" NX5 is resolved by IRAC (see Fig.fTsl. 

^ NXIO is located at the tip a prominent pillar structure in the M17 H II region. IJ02I resolve this source 
into two YSOs of comparable brightness in the near-IR, a Class I and a Class II/III candidate. The Class 
I source likely dominates the mid-lR emission. 



lifetimes. 

Assuming a steady rate of star formation, the rela- 
tive number of Stage II versus Stage O/I YSOs gives the 
relative lifetime of each Stage. But we must take care 
to compare populations with similar mass distributions. 
IVIost YSOs fit by Stage I models have steep, Class I 
spectral indices, and even low-mass Stage O/I YSOs can 
be bright in the mid-IR. Consequently, Stage O/I sources 
disproportionately populate the YMF for Af* < 3 IVIq 
(Fig.[IT]), with the 12 Stage O/I candidates in the YJVIF 
(Table [2]) overrepresented at 1 ]VIq< < 2 JVIg com- 
pared to Stage II candidates. For > 3 IVIq, the YJVIF 
contains 5.6 times as many Stage II as Stage O/I candi- 
dates, while this ratio is 2.6 for the entire mass range. 
The median Stage II disk lifetime in our sample is there- 
fore tu = 5.6tA ~ 0.4 IVIyr. The above calculations ex- 
clude the single Stage III candidate and 18 candidate 
YSOs with ambiguous Stage determinations in Table [21 
Nearly all of the ambiguous candidates have model fits 



that are divided between Stage O/I and Stage II mod- 
els. In the limiting case where all ambiguous sources are 
treated as Stage II candidates, we have to — 8.2i^ < 0.6 
JVIyr. Intermediate-mass YSOs therefore lose their disks 
in less than half of the corresponding '--^l JVIyr timescale 
observed for their solar-mass counterparts. This age dif- 
ferential tends to steepen the slope of the YMF com- 
pared to a disk-unbiased IJVIF, hence the intermediate- 
mass YMF probably underestimates the true number of 
low-mass YSOs. 

The characteristic age of the source population in the 
YJVIF is tA -\- to ^ 0.5 IVIyr. Using the stellar mass esti- 
mated from the IMF model, we arrive at a star formation 
rate (SFR) of 0.0016 Mg yr'^ in the extended target 
field, which applies to M17 EB and the extended outer 
regions of the M17 molecular cloud. Due to the poten- 
tial incompleteness and the bias against diskless sources 
inherent in the intermediate-mass YMF, this value must 
be regarded as a lower limit. 
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Fig. 12.— Image showing part of the interface between M17 EB and MC G15.9-0.7 (IRAC [4.5]=blue, IRAC [8.0]=green, 
MIPS [24]=red). The cyan contours are '^^CO emission at u = 19 km s~^. The rim of the bubble M17 EB is delineated by the 
bright green filament of 8 /xm emission crossing from the lower-left to the upper-right corner of the image. Candidate YSOs are 
overplotted as in Fig. 1101 Two candidate OB stars associated with diffuse 24 ^m emission are marked by white diamonds and 
numbered as in Tabled The candidate Stage O/I YSO G015. 6653-00.4989 (E55 in Table [I} is highlighted, and the model fits 
to its SED are plotted. The [4.5] fiux was used as an upper limit to the fitting and is plotted as a triangle. The dashed curve 
shows the reddened photosphere of the central star as it would appear in the absence of circumstellar material. 



4.2.3. Triggered Star Formation m MC G15.9-0.7? 

A concentration of star formation activity is observed 
in the vicinity of an IR dark cloud (IRDC) complex cen- 
tered at {l,b) = (15.55°, -0.45°). This region is shown 
in a composite GLIMPSE and MIPSGAL image in Fig- 
ure 1121 The bright, filamentary 8 /im emission slicing 
diagonally through this image, part of which makes up 
a previously identified optical reflection nebula, outlines 
the rim of M17 EB where it intersects MC G15.9-0.7 (Fig. 
[1]). Because significant concentrations of candidate YSOs 



are observed on the bubble rim while no group of candi- 
date YSOs is reliably detected in the adjacent regions of 
MC G15.9-0.7 (Fig.[Tni), the situation is strongly sugges- 
tive of star formation triggered by the expansion of M17 
EB. The morphology of the IRDCs, bright 8 /im emis- 
sion, and CO emission provides evidence for increased 
gas density resulting from the compression of the cloud. 

The SED of G015. 6653-00.4989, a bright 24 /xm source 
in MC G15.9-0.7 (Fig. [121), is very similar to the SEDs 
of several massive pr otostars associated wit h molecular 
outflows identified by I Shepherd et all (|2007f ). It is char- 
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acterized by excess emission at 4.5 /im (an outflow sig- 
nature) and a very red IR spectral index peaking near 
70 /im (this object is the brightest 70 /im point source 
detected in the M17 target field). The models fit to this 
source (E55 in Table [2]) indicate a luminous candidate 
protostar with mass 7 Mq < Mi, < 8 Mq and luminos- 
ity 400 L0 < Ltot < 10'^ L0 accreting from its cir- 
cumstellar envelope at a high rate of 2 x 10"'' Mq yr"'^ 
< Mcnv < 8 X 10"^ Mq yr~^. The accretion age of this 
candidate protostar is tA < 0.04 Myr. A nearby resolved, 
compact, red IR source (in the white box in Fig. [12]) is 
associated with a second 70 iim point source. The SED 
of this source cannot be fit well by any YSO model, and 
we speculate that it may be either a luminous binary 
YSO or a compact cluster. Both sources are apparently 
located in the cloud region compressed by M17 EB. 

Also shown in Figure [12] are 2 compact, resolved 
sources of 24 /im emission (30"~60" in size) that are not 
associated with any candidate YSO. This emission sug- 
gests a localized region of heated dust, and the most 
likely heat source is a massive young star or cluster. 
Five such 24 fim sources are observed in the Ml 7 tar- 
get field, and none are found in the control field. At the 
center of every 24 fj,m source lies a GLIMPSE Archive 
source with an SED consistent with a hot stellar photo- 
sphere . Following the method pioneered bv lWatson et al.l 
([2008( 1 ■ we estimated the spectral types of the se stars 
Iw employing the RW07 model fitting tool to fit iKurucd 
([1993) stellar atmospheres (including hot photospheres 
up to Toff = 50, 000 K) to available BVJHKs+IRAC 
photometry. We found that all of the stars associated 
with compact 24 //m emission have luminosities consis- 
tent with B or late O stars, assuming they are associ- 
ated with M17 at either a distance of 1.6 kpc or 2.1 
kpc (see Table [5]) . Because these stars exhibit little or 
no IR excess emission in the 2MASS/GLIMPSE bands, 
they either have already cleared out the inner regions of 
their circumstellar disks or have dissipated their disks 
entirely. The observed 24 /zm emission could arise from 
remnant disk material or heated dust in the surrounding 
molecular cloud at distances of ~2-5 x 10^ AU from the 
central star. This interpretation suggests that some of 
these stars may ionize low-luminosity compact H II re- 
gions that_Jiave_jiot_been det^^ by previous surveys 
(e.g. iWood fc Churchwel]|[l989h . None of these sources 
is det ected by the high-resolutioii COR NISH survey at 
6 cm ([Purcell. Hoare. fc Diamondll2008D . indicating that 
any gas ionized by these B stars is diffuse enough to be re- 
solved out by the VLA in the B configuration. This rules 
out the presence of hypercompact or ultracompact H II 
regions around the stars. The most extended of these 
24 /xm sources, associated with star 24 fFi g. [T2l Table 
5]), is associated with 11 cm radio emission ([Reich et al.l 
I990| ) and enclosed by a partial ring of 8 ^m emission at 
a projected distance of 0.5 pc from the star, suggesting 
a small H II region and PDR. 

Whatever the interpretation of the physical origin of 
the 24 /Ltm sources, the circumstellar environments of the 
associated OB stars are significantly more evolved than 
those of the candidate YSOs, which exhibit IR excess 
emission at wavelengths <10 /zm. The comparison of 
star 25 with YSO E55, the candidate protostar in Figure 
1121 is particularly interesting. These two sources have 



comparable masses (~10 M©; E55 is still gaining mass 
through accretion) and are apparently located within the 
same molecular cloud, but they occupy strikingly differ- 
ent phases of early massive stellar evolution, suggesting 
a difference in age comparable to the disk destruction 
timescale derived above, tD 0.4 Myr. This difference 
in age could be explained naturally by triggered star for- 
mation. The more evolved B star lies in the midst of the 8 
/xm emission defining the inner rim of M17 EB, while the 
young candidate protostar is located in a region of the 
molecular cloud that has only recently been overtaken 
by the expanding bubble. The expansion velocity of the 
bubble implied by this scenario is slow, a few km s~^ at 
most, which is consistent with the non-detection of an 
expansion signature in the CO line (§4 above). 

4.3. The M17 Molecular Cloud and H II Region 

The contours in the top panel of Figure [8| show a 90- 
cm V ery Large Array (VLA) radio image ([Brogan et al.l 
l2006f ) delineating the extent of the ionized gas associ- 
ated with the main M17 H II region. The densest part 
of the M17 molecular cloud is located in M17 South (see 
Fig. [1]) along the interface between the H II region and 
the bright '^■^CO (J = 2 ^ 1) emission (colored red- 
white in Fig. [8]). The total gas mass (molecular and 
atomic) contained in the main M17 molecular cloud is 
^6 X lO'' M0, found by integrating the column density 
traced by the CO emission. The H II region occupies a 
cavity cleared of molecular gas between M17 North and 
M17 South. The morphology of the CO emission from 
M 17 closely resembles the maps o f the region obtained 
bv [Wilson. Howe, fc Balogh[ ([1999[ ). but the older maps 
have slightly lower resolution and do not include M17 
EB. The cavity in the molecular cloud has been sculpted 
by the winds and radiation of the O stars in NGC 6618 
(PSC07). 

4.3.1. Prominent Massive Stars 

The spectral classifications derived from SED fits to 
27 prominent OB stars in t he M17 target fleld are pre- 
sented in Table [3 We fit Kurucd (ji993) stellar atmo- 
spheres to the broadband fluxes of 11 of the 16 known 
O stars in NGC 6618 that have reliable, unsaturated 
fluxes in the GLIMPSE Archive to estimate their spec- 
tral types. Although all known O stars in NGC 6618 
except No. 8 have previously reported spectral types 
from optical and near- I R spectroscopy and photometry 
(lOgura fc Ishidal [197^ [Chmi. Elsasser. fc Neckell [19801 : 
[Hanson. Howarth. fc Conti| [1997: H0"8l). mid-IR photom- 
etry offers the advantage of a relatively extinction- 
independent determination of the luminosity in the 
Rayleigh- Jeans regime of the stel lar SEP. We incorpo- 
rated available B, V photometry ([Ogura fc Ishida[ri976[ : 
[H08f ) into each SED to obtain more reliable estimates of 
Ay toward each star than would be possible from the 
2MASS photometry alone. Because the optical-IR SED 
shape does not vary significantly for stars earlier than B3, 
the Ay derived from our SED fltting is independent of 
adopted distance. As shown in Table [5] our SED-based 
luminosity classification at 1.6 kpc agrees with previously 
reported spectral types to within 2 subclasses for 8 of the 
O stars analyzed in NGC 6618. This might seem surpris- 
ing, given that the most recent distance determination 
found that M17 cannot be significantly nearer than 2 
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TABLE 5 

Selected OB Stars Associated with M17 



Index 




I 


b 


Cataloged 


Spectral Type from SED 


Av 


No.'' 


Name'' 


(deg) 


(deg) 


Sp. type': 


d = 1.6 kpc 


d = 2.1 kpc 


(mag) 


O Stars in NGC 6618 


1,2=^ 


GEN la,b 


15.0562 


-0.6884 


04+04 V 








3 


GEN 43 


15.0533 


-0.7045 


03-05 V 


04-05 V 


05-f05 V 


12.3 


4 


GEN 2 


15.0731 


-0.7004 


05 V 


05.5 V 


05.5-I-07 V 


5.2 


5 


GEN 37 


15.0559 


-0.6883 


03-06 V 








6 


OI 345 


15.0110 


-0.7020 


06 V 








7 


GEN 18 


15.0812 


-0.6569 


06-08 V 


06 V 


06+07.5 V 


7.6 


8 


M17-S3 


15.1032 


-0.6487 




07 V 


05 V 


11.2 


9 


GEN 25 


15.0673 


-0.6879 


07-08 V 


09 V 


06.5 V 


8.0 


10 


OI 352 


14.9945 


-0.7486 


08 V 


05.5-f05.5 V 


multiple": 


~7 


11 


OI 174 


15.1325 


-0.5257 


09 V 


06 V 


04 V or 06+07 V 


7.0 


12 


GEN 3 


15.0658 


-0.7084 


09 V 


07 V 


05 V or 07+07 V 


3.7 


13 


GEN 16 


15.0748 


-0.6460 


09-B2 V 


O9.5-B0.5 V 


08 V or BO-I-BO V 


5.9 


14 


GEN 61 


15.0594 


-0.6884 


09-B2 V 


09 V 


06.5 V or 09.5-I-09.5 V 


9.9 


15 


GEN 27 


15.0433 


-0.6950 


09 V 


Bl V 


09 V 


~9.5 


16 


GEN 31 


15.0759 


-0.6534 


09.5 V 








Candidate OB Stars in NGC 6618PG 


17 


BD-16 4831 


15.3213 


-0.7758 


0+ 


05.5 V or 09 III 


04 V or 07 III 


4.4 


18« 


BD-16 4826 


15.2604 


-0.7263 


05 


05.5 V or 09 III 


05 V or 06.5 III 


3.9 


19 


BD-16 4822 


15.2244 


-0.6796 


B... 


08.5 V 


06.5 V 


2.9 


20 


TYC 6265-1174-1 


15.1882 


-0.7598 


BO 


Bl V 


09 V 


3 


21 


HD 168585 


15.2855 


-0.7498 


B7-B8 II 


Bl V 


09 V 


1.6 


Candidate OB Stars Associated with 24 /xm Emission in M17 EB 


22S 


ISOGAL YSO 


15.3320 


-0.7174 




late B V 


B1.5 V 


2.1 


23^ 


TYC 6265-347-1 


15.3583 


-0.6543 


O... 


B2 V 


Bl V 


2.6 


24 


BD-15 4928 


15.4801 


-0.3889 


B 


B0.5 V 


08 V 


2.8 


25 




15.5702 


-0.4895 




Bl V 


B0.5 V 


~4.8 


26 


TYC 6265-2079-1 


15.5900 


-0.8268 


O... 


B4 V 


B2 V 


1.7 



Note. — All O stars in NGC 6618 are Chandra point sources (BFT07). The only star in NGC 6618PG that has been observed by 
Chandra, No. 20, is also an X-ray point source. 

^ OB stars in each of the first two groups are listed in order of decreasing luminosity, while the candidate B stars in the final group are 
listed in order of increasing Galactic longitude. 

^ CEN is lChini. Els asser. Sz Ncckel' (fT980|), and OI is lOgura & I shida (197^. SWB M17-3 is a recently-discovered O star driving a stellar 
wind bow shock ijPo vich et al. 2008). Other identifiers are from SIMBAD. 

^ Spectral types for NGC 6618 are from OI, CEN. lHanson. Howarth, & Contil ( I1997D . and lHOSt others are from SIMBAD. 

CEN la,b refer to the main components of an apparent 04-I-04 visual binary that are saturated and marginally resolved by IRAG. 
GEN la,b each are spectroscopic binary systems with nearly equal-luminosity components ( |H08| ). 

^ At 2.1 kpc, the SED fits to OI 352 return a radius larger than that of an O supergiant. It may be an unresolved multiple O system 
(see text). 

^ BD-16 4826 is a ROSAT point source jDunne et 'an[2003l) . 

^ This source, ISOGAL-P J182108. 8-155658 from the catalog of IFelli et aP ||2002|1 . exhibits a marginal IR excess at 5.8 fim and is not 

detected as a point source at 8.0 fim. 

^ This source exhibits an IR excess at 8 /j.ni. 



kpc (|H08l ). However, previous photometric determina- 
tions of spectral type assumed that each source was a 
single star, and our method shares the same systematic 
bias with respect to distance. Even spectroscopic de- 
terminations have generally neglected multiplicity, and 
the resulting spectrophotometric distance estimates con- 
tributed to the wide adoption of the 1.6 kpc distance. 

Revising the distance from 1.6 to 2.1 kpc multiplies 
the luminosity of every star in Table [5] by 1.7. This 
factor of ~2 is suggestive. The majority of O stars in 
clusters have been observed to b e in near- cqual-ma ss bi- 
nary systems (Sana ct al. 2008; Mason et al. 2009). IH08I 
found that each component of the historical 04-1-04 bi- 
nary CEN 1, the most luminous star system in M17, is 
itself a spectroscopic binary with approximately equal- 
luminosity components and suggested that several other 
known O stars in NGC 6618 could also be unresolved 
binaries. Assuming the stars in Table [5] are at 2.1 kpc, 
we find that the majority of O stars in NGC 6618 are 
far more luminous in the mid-IR than expected from the 
cataloged main-sequence spectral types. Three of these 



stars, CEN 43, CEN 2, and CEN 18, each have mid-IR 
luminosities exceeding an 03 V star. Three others, 01 
174, CEN 3, and CEN 61, are unhkcly to be single main- 
sequence stars because their SED-derived spectral types 
are > 2 subclasses earlier than their cataloged spectral 
types (Table [5]). Because NGC 6618 is very young and 
even 04 V stars take ^4 M yr to evolve off the main se- 
quence (|Bressan et al.|[T993f) . it is unlikely that the mas- 
sive stars in NGC 6618 are giants. We therefore give bi- 
nary classifications for the over-luminous stars in Table 
[HI with spectral types estimated for each component as- 
suming (1) that the components have nearly equal masses 
and (2) that their combined light yields the cataloged or 
SED-derived spectral type for a single star at 1.6 kpc. 
We stress that our spectral typing for each component of 
the binary systems is not definitive, but the high fraction 
of candidate binary systems is a robust result. 

At 2.1 kpc, one source, 01 352, has mid-IR luminosity 
equivalent to a binary 03-1-03 V system. If this were the 
true classification, then 01 352 would compete with CEN 
la,b for domination of the H II region, but in fact 01 352 
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is located in an unremarkable part of the nebula, outside 
thecentral wind-dominated cavity (jTownslev et al.l[2003t 
IPSCOTI ). We therefore speculate that this star, like CEN 
1, could be a multiple system of 3 or more later- type 
0-dwarfs. BFT07 report an unusually hard (4.4 keV) 
thermal component from their plasma fits to the X-ray 
spectrum of 01 352. Strong shocks produced by coUiding- 
wind binaries are the most promising explanation for O 
stars with >1 keV hard X-ray emission. CEN la,b, CEN 
43, CEN 2, and CEN 18 aU produce hard X-ray emission 
(BFT07), and all are identified as candidate binaries by 
H08 and/or by our SED fitting, with primaries of type 06 

V or earlier. CEN 25, a solitary 06.5 V star, lacks hard 
X-ray emission. In contrast, 01 174, CEN 3, and CEN 
61 are likely binary systems with primaries of type 06 

V or later, and these systems have only soft-component 
X-ray spectra. Evidentally hard X-ray emission is not a 
ubiquitous property of massive binary systems, but this 
phenomenon may be more common in earlier-type bina- 
ries. CEN 16 seems anomalous, for although it boasts 
a very hard X-ray spectrum (3.5 keV), it is among the 
least luminous candidate O stars, equivalent to an early 
B binary system (Table [5|) . 

4.3.2. Comparison of X-ray Sources with GLIMPSE YSOs 

The white box in the bottom panel of FigurefTOloutlines 
the field observed to date by Chandra (Townsley et al. in 
prep., Broos et al. in prep.). BFT07 presented 886 X-ray 
point sources extracted from the first Chandra Advanced 
CCD Imaging Spectrometer (ACIS) observation of M17. 
By scaling the X-ray luminosity function (XLF) of the 
Orion Trapezium Cluster to match the XLF of the M17 
sources, BFT07 estimated a total young stellar popula- 
tion of 8,000-10,000 members (to a limiting mass of 0.1 
Mq) in the NGC 6618 cluster and its immediate vicinity 
in M17. BFT07 assumed a distance of 1.6 kpc, and scal- 
ing their XLF to 2.1 kpc increases the total population 
to 11,000-14,000 stars and YSOs. 

Of the 195 candidate YSOs in our M17 subsample (Ta- 
ble [T]), only 34 fall within the 17' x 17' field of the ACIS 
observations. This number would be many times higher 
if not for the presence of extremely bright diffuse mid-IR 
emission from the H II region and PDR that drastically 
reduces the point-source sensitivity in all IRAC bands. 
The source density in the GLIMPSE Archive within the 
ACIS field analyzed by BFT07 is 107,000 deg-^. This is 
only ~50% of the average density of Archive sources in 
the larger M17 target field (see Table [T]), despite the fact 
that the ACIS field contains the dense center of the NGC 
6618 cluster. Because of this incompleteness, we cannot 
attempt any meaningful quantification of the YSO pop- 
ulation most directly associated with NGC 6618. X-ray 
observations, by contrast, do not suffer from the dele- 
terious effects of bright backgrounds, and moreover the 
(near-axis) spatial resolution of ACIS is higher than that 
of IRAC. Because very different physical mechanisms 
are responsible for X-ray versus IR emission from young 
stars, we are interested in learning whether the IR-excess 
sources selected for our YSO sample represent a subsam- 
ple of the X-ray-selected population or a distinct popu- 
l ation. 

iPreibisch et al.l ()2005l ) performed a detailed analysis 
of the X-ray source population observed as part of the 
Chandra Orion Ultradeep Project (COUP). They found 



that the 838.1-ks COUP observation detected 100% of 
the known young F- and G-type young stars and YSOs 
in Ori on but only ~50% of A-type stars. IPreibisch et al.l 
(|2005f) concluded that the X -ray emission observed in 
the low-mass stars in the COUP sample is consistent 
with an origin in stellar coronae analogous to the so- 
lar corona, but at a considerably higher level of activity. 
The disappearance of solar-type convective envelopes in 
more massive stars on the main sequence could explain 
the drop in t he COUP comp leteness for intermediate- 
mass stars. IPreibisch et al.l (|2005,) also observed that 
X-ray emission in T Tauri stars is partially suppressed 
by accretion. As we have shown, our IR-selected sample 
of candidate YSOs is highly sensitive to intermediate- 
to high-mass objects that are still accreting via circum- 
stellar disks and/or envelopes. We therefore expect that 
our YSO sample complements the (much larger) X-ray- 
selected sample, as we preferentially select the YSOs that 
are the most likely to be missed by Chandra. 

We have correlated our final YSO sample with an up- 
dated X-ray source catalog covering the entire area ob- 
served to date by Chandra/ ACIS, with 400 ks total inte- 
gration time (Broos et al. in prep). Out of the 34 candi- 
date YSOs in the overlap region, 14 have X-ray counter- 
parts. The 95% confidence intervals for M^, and Ltot, 
along with the evolutionary Stage derived from the well- 
fit (Equation [T]) YSO models, are given in Table ID for 
each source. The sources in Table |3I exhibit no obvious 
trends in mass, luminosity, or evolutionary Stage with 
which to distinguish candidate YSOs with X-ray coun- 
terparts from those without. The majority of the candi- 
date YSOs in Table|4]are fit with models of intermediate 
to high stellar mass corresponding to spectral types A 
and B on the main sequence. There may be a slight 
preference for less-evolved sources with higher accretion 
rates in the X-ray quiet population in Table HI but this 
is not significant, given the small sample size. Intrinsic 
X-ray variability may play a role in determining which 
YSOs are detected. It is also possible that some of the 
luminous YSOs with ACIS detections are themselves X- 
ray quiet but possess binary companions that are X-ray 
bright. There is no evidence that the 34 candidate YSOs 
listed in Table^ represent a distinct population from the 
X-ray-selected sources. 

Two interesting candidate massive YSOs are located 
near the peak molecular gas density of Ml 7 South and 
the interface with the H II region, shown in Figure [131 
One, the KW Object, is a binary system. A candidate 
Herbig Be star is the more lumi nous component (KW-1 ) 
and dominates the IR emission djOl : IChini eral]f2004ah . 
The lower-luminosity component (KW-2) is X-ray bright 
(BFT07) but unresolved in our IR survey data, hence the 
KW Object is considered an X-ray detected source (X4) 
in Table H KW-1 is saturated in the GLIMPSE and 
MIPSGAL images, so models were fit to its 2MASS and 
MSX fluxes only. The model fits indicate a central star 
of ~ 10 Mq and Toff ~ 26, 000 K, consistent with 
an earl y B star an d in go od agreement with the find- 
ings of IChini et al.l (|2004al) . The models indicate that 
the star may have a substantial circumstellar disk (0.01 

Microshocks in OB stellar winds are believed to be the primary 
source of soft X-ray emission from massive stars; this very different 
from the emission mechanism operating in low-mass stars. 
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Fig. 13. — Region containing the peak molecular gas density 
in M17, near the H II region-PDR interface. The image com- 
bines IRAC [4.5] (blue) and [5.8] (green) with MSX E (21.3 
^m; red). The cyan contours are the HHT ^^CO moment 
map from the top panel of Fig. [S] Candidate YSOs are 
overplotted as circles using the color scheme of Fig. [S] The 
uhracompact H II region M17 UC-1 is marked. The KW 
Object (left) and neighboring resolved YSO [right] are also 
highlighted, and the model fits to their SEDs are plotted as in 
Fig. 1121 Fluxes used as upper or lower limits for the model fit- 
ting are plotted as downward- and upward-pointing triangles, 
respectively. 



Mo< A/disk < 0.3 M0) with a large inner hole of ra- 
dius ~40 AU (^10 times the dust sublimation radius). 
The inner hole is responsible for the precipitous rise in 
the near- through mid-IR SED (Fig. [0)1 and indicates 
that the central star has ceased accreting material from 
the disk. Our model SEDs produce a broad 9.7 /xm sili- 
cat e feature in emission, but ISOCAM spectra presented 
by IChini et al.l ( 2004a) show that the "bumps" in the 
SED near 8 and 11 /^m are actually due to PAH emis- 
sion features excited by the hot, central star. We have 
therefore used the MSX A and C bands as upper lim- 
its to the SED fitting, since these bands contain PAH 
features that increase their flux with respect to the con- 
tinuum provided b y the YSO rn odels. The 1-D envelope 
model employed bv lChini et al.1 (|2004al predicts the sili- 
cate feature in absorption and does not explicitly include 
a disk. These authors report a total circumstellar mass 
of 10 Mq within 0.1 pc from the star. The total mass 
and envelope radius are very difficult to constrain in the 
absence of mm/sub-mm data, hence we cannot confirm 
these reported values. 

The other prominent candidate YSO highlighted in 
Figureim G014.9791-00.6651, is resolved by IRAC. This 
source, NX5 in Table [H is not detected in X-rays. Like 
KW-1, NX5 is massive [AU ~ 10 Mq), but it may be 
somewhat less evolved. The parameters of the 47 well-fit 



models poorly constrain the circumstellar environment 
of NX5, allowing for a wide range of both Mdisk and 
A/cnv, so the evolutionary Stage is ambiguous, split be- 
tween 60% Stage II and 40% Stage O/I models. The 
resolved outer envelope of this star is large, with radius 
'^10^ AU, and the SED exhibits the signatures of a PAH 
spectrum. As with our modeling of KW-1, we set the 
bands most affected by PAHs as upper limits for the 
model fitting, and the majority of well-fit models clus- 
ter around ToH' ~ 23,000 K, hot enough for the star to 
produce sufficient UV radiation to excite PAHs in its 
circumstellar material. Together with M17 UC-1, the 
near by hypercompact H II region that is excited by a BO 
star (|Sewilo et al.ll2004D . KW-1 and NX5 are 3 very lumi- 
nous YSOs with very similar stellar masses but with very 
different circumstellar environments. This demonstrates 
that, in spite of their close proximity to each other, these 
3 sources are not precisely coeval. 

4.3.3. Challenges of Identifying Cluster Members and 
Determining the Age of NGC 6618 

The fraction of stars exhibiting IR excess emission due 
to circumstel lar disks is a diagnostic of the age of embed- 
ded clusters ()Haisch. Lada. fc Ladal[2001b[ ). BFT07 find 
that only ^^10% of the sources in their X-ray selected 
sample exhibit IR excesses. This fraction is very low 
in comparison to recent results using IR-selected source 
samples. |J02| found a JHK excess fraction of 53 % in 
the central regions of the NGC 6618 cluster, and IHOSI 
reported a 74% excess fraction from deeper JHKL ob- 
servations of a similar field. Some variation in excess 
fraction is inevitable due to differences in sample selec- 
tion, and excess fractions measured in the JHKL color 
plane are expected to be higher than those in the JHK 
color plane. The low excess fraction reported by BFT07 
is partially explained by selection biases against embed- 
ded sources and against low-mass sources, both of which 
are fainter in X-rays. To gain depth beyond the 2MASS 
and GLIMPSE catalogs, BFT07 correlated their X-ray 
sample with SIRIUS JHKs photometry presented by 
J02, hence the difference between the excess fractions re- 
ported by these two works reflects the different selection 
criteria and analysis techniques rather than differences 
in IR photometry. 

The large discrepancy between reported IR excess frac- 
tions reflects a more basic issue: The M17 H II region is 
a difficult place to search for YSOs. The diffuse IR emis- 
sion is bright at all wavelengths with a steep red spec- 
tral index and exhibits complex morphological structure 
on multiple spatial scales (|PSC07D . This makes reliable 
mid-IR photometry difficult even for bright sources and 
impossible for faint ones. Nebular structures can be mis- 
taken for IR excess emission associated with stars. As 
iPovich et al.l (|2008( ) note, two of th e high- mass Class I 
sources identified by iNielbock et al.l (|2001h are actually 
mid-IR emission from stellar-wind bo w shocks tha t could 
be mistaken for circumstellar disks. iLada et "ahl ()1991l ) 
claimed that the majority of high-mass stars in NGC 
6618 display JHK excess, but this was disputed bv lHOSl . 
Among the O stars in NGC 6618 detected by GLIMPSE, 
all are detected at wavelengths as long as 4.5 /im and 
many are detected at 8.0 /im, but we can state with 
high confidence that none displays an IR excess, with 
the possible exception of star No. 8 in Table [5l which is 
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assoc iated with a third mid-IR bow shock (jPovich et al.l 
[2OOI . 

Near-IR photometry can be more rehable, because 
the diffuse emission is less bright. But JHKL obser- 
vations are chaUenged by the large spatial variations 
in interstellar extinction across the M17 cloud (see the 
Ay values in Table [51 for example). While the ex- 
tinction due to molecular gas in the M17 North and 
South cloud components is generally Ay > 15 mag, 
and the dense core of M17 South produces Ay > 90 
mag, the average extinction through the central cavity 
is Ay < 15 mag (Fig. [S]). The average extinction ob- 
served toward the ionizing stars of t he M17 H II region 
is Ay ^8 mag (Table [51 Chini, Elsa sser. fc NeckellllQSOl : 
iHanson. Howarth. fc Contiill9971 : IHOl ). The fore ground 
extinction to the cloud is Ay ~ 2 mag (|Pellegr ini et al.l 
I2007HH08D . We thus deduce that the column density of 
the M17 cloud is divided into approximately equal halves 
behind {Ay 7 mag) and in front {Ay ~ 6 mag) of the 
H II region. Including foreground extinction, the total 
extinction along a typical line-of-sight through the cav- 
ity is Ay < 17 mag, or Ak < 2 mag. Hence background 
stars are detectable even in relatively shallow near- and 
mid-IR observations along the line of sight through the 
central cavity and the ionizing cluster. In contrast, ad- 
jacent dense regions of the M17 North and South cloud 
components produce sufficient extinction to screen out 
background stars against even deep near-IR observations. 
Because of this complex extinction morphology, there is 
no comparable nearby field that can be used to estimate 
the background contamination in observations of NGC 
6618. Any attempt to measure accurately the disk frac- 
tion or low-mass IMF in NGC 6618 must overcome this 
significant challenge. 

A major strength of X-ray observations is the ability to 
separate young cluster members, which are 2-3 orders of 
magnitude brighter in X-rays than older main-sequence 
or evolved stars, from the field star population (BFT07 
and references therein). The contamination from unas- 
sociated sources in the X-ray sample of the M17 young 
stell ar population is therefore very low, <5% (BFT07). 
IJO2I recognized the difficulty in using their off-cluster field 
to correct for source contamination and instead used a 
Galactic disk population synthesis model. While we are 
not in a position to verify their results, the conservative 
upper limit of 3 Myr on the cluster age derived by J02 is 
reasonable. H08 assumed that all stars with IR excess, 
and by extension all sources with reddening equivalent 
to Ay > 10 mag, were cluster members. H08 therefore 
concluded that the M17 molecular cloud behind NGC 
6618 was opaque. Both the assumption and the conclu- 
sion were not valid. The extinction through the central 
cavity is relatively low {Ak < 2 mag) and the deep obser- 
vations used by H08 detected stars as faint as K — 20.1 
mag and L = 15.6 mag. The sample of H08 therefore 
contains a potentially large fraction of background stars, 
casting doubt upon the reported disk fractions and the 
correspondingly young age of 0.5 Myr claimed for the 
M17 stellar population. 

It is tempting to conclude that the inclusion of disk- 
less field stars in a magnitude-limited sample of the 
M17 young stellar population would decrease the mea- 
sured IR excess fraction. But stars with disks are 
not well-separated from reddened stellar photospheres in 



the JHK or even in t he JHKL color-color diagrams 
(jWhitnev et all I2003U iRobitaille et"all l2006l ). Color- 
color spaces or SEDs that include photometry at 4.5 fim 
and longer wavelengths provide a significantly more re- 
liable means for distinguishing circumstellar extinction 
from interstellar extinction. Of the 64 candidate X- 
ray emitting protostars identified by BFT07 (their Ta- 
ble 7) most of them on the basis of the K — [3.6] color 
(analogous to K — L), only 3 were selected as YSOs by 
our more conservative multiband SED fitting (see Ta- 
ble [4]), but 16 were well-fit by reddened stellar photo- 
spheres (the remaining sources are detected in too few 
2MASS/GLIMPSE bands to be fit rehably with SED 
models). Highly-reddened contaminating sources that 
are intrinsically red, such as the numerous red giants 
in the Galactic plane behind the M17 H II region, are 
the most difficult to separate from YSOs on the basis 
of JHKL colors alone. If correction for contamination 
is not done, then highly-reddened, background red gi- 
ants may provide a significant number of false IR excess 
sources. This would artificially increase the derived disk 
fraction, leading to an underestimate of the cluster age. 

We are left without a strong constraint on the age of 
the M17 young stellar population. There is a population 
of PMS stars lacking inner disks (BFT07), and there is 
also ample evidence for ongoing star formation (Nielbock 
et al. 2001, J02, this work). The preponderance of obser- 
vational evidence indicates that M17 has not formed all 
of its stars simultaneously, but has experienced continu- 
ous or episodic star formation over a period of <3 Myr. 
Sources like the candidate massive YSOs and M17 UC-1 
shown in Figure [T3l are almost certainly younger than 0.5 
Myr, the average age of the extended YSO population 
(§4.2 above). To minimize the impact of the lingering 
age uncertainty, we adopt a mean age of ~1 Myr for the 
young stars and YSOs in NGC 6618 and the surrounding 
molecular cloud. 

We can use the BFT07 results to make an estimate 
of the star formation rate (SFR) that has produced the 
ionizing cluster of M17 and its entourage of YSOs. Scal- 
ing t he stellar initia l mass function (IMF) measured by 
Muen clTet al.l ([.2002) for the Orion Trapezium Cluster to 
the 11,000-14,000 stars in M17, we arrive at a total mass 
in stars of 8,000-10,000 M©. Stars have been forming in 
M17 at an average rate of 0.008-0.01 Mq yr'^ 

5. A HISTORY OF PROPAGATING MASSIVE STAR 
FORMATION 

We have shown that the major sites of star formation 
in the extended M17 region are located around the rim 
of the large bubble M17 EB. We will now examine the 
structure and possible origins of Ml 7 EB and show that 
the beginnings of a plausible star formation history of 
M17 are found within this bubble. 

5.1. The Origin of Mil EB 

The most likely production mechanisms for a bubble 
the size of M17 EB are the expansion of an H II region 
or a supernova explosion, both of which are expected 
to occur in regions of massive star formation. We fa- 
vor the H II region interpretation. The morphology of 
the mid-IR emission outlining M17 EB is similar to that 
of oth er 8 /im bubbles outlining the PDRs of H II re- 
gions (|Churchwell et al.ll2006H Watson et al.ll2008f) . One 
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of the best ways to discriminate between H II regions and 
supernova remnants (SNRs) is by the radio continuum 
spectral index. H II regions are filled with thermal emis- 
sion from ionized gas and exhibit a flat spectral index 
for optically thin radio emission, while SNRs arc often 
characterized by steep, non- thermal spectral indi ces from 
synchrotron emission (e.g. iBrogan et all l2006f) . M17 
EB is included in two single-dish radio Galactic plane 
surveys of comparable spatial resolution. Images at 3 
cm from the Nob eyama 45-m telescope (HPBW=2.7'; 
iHanda et allllOSTl ) and at 11 cm from the Effe lsberg 100- 
m telescope (HPBW=4.3'; iReich et al.lll990t )^^ show a 
plateau of radio continuum emission that fills the half 
of M17 EB farthest from M17 (other half is dominated 
by emission from the M17 H II region itself) and drops 
sharply at the location of the 8 fim emission outlining 
M17 EB. Comparison of the 3 cm and 11 cm surface 
brightnesses at various points in the interior of M17 EB 
shows that the radio continuum flux density is approxi- 
mately constant (with a typical value of ~50 mJy) over 
this wavelength range, consistent with a flat, thermal 
emission spectral index. 

One caveat regarding the interpretation of the radio 
emission within M17 EB should be kept in mind. The 
non-thermal emission from SNRs eventually fades (on 
the order of a few 10^ yr) to the point of becoming indis- 
tinguishable from the background. The residual thermal 
emission from the parent H II region lasts much longer. 
Hence, we cannot rule out the possibility that a super- 
nova has occurred within M17 EB. But because the ob- 
served thermal radio emission completely fills the bubble, 
there is no need to invoke a supernova explosion to ex- 
plain its origin. 

A search of SIMBAD objects with coordinates inside 
M17 EB reveals a group of 5 stars with OB spectral clas- 
sifications. Although a few of the se stars have proper 
motion estimates (Hog et aLll2000l ). none have parallax 
distances. The spectral types of these sources listed in 
the SIMBAD database are of unknown provenance and 
questionable quality, so while we list the catalogued spec- 
tral types for reference in the lower 2 sections of Table [HI 
we trust the classifications only to the extent that the y 
identify hot stars. Again fol lowing I Watson et al.l (|2008f ). 
we fit SEDs generated from iKurucd (|1993l ) stellar atmo- 
spheres to available B, V, 2MASS, and GLIMPSE fluxes 
of these stars, and we find that their IR luminosities are 
consistent with very early spectral types at either 1.6 kpc 
or 2.1 kpc (Table [5|). If the stars were at a significantly 
different distance, their SEDs would be inconsistent with 
early-type stars. We therefore suggest that these stars 
are the 5 most massive members of NGC 6618PG, the 
ionizing cluster driving M17 EB f Figs . [HI and [T5)l . 

The OB stars of NGC 6618PG are significantly less 
reddened {Ay < 5 mag) than the O stars in NGC 6618, 
as expected since NGC 6618PG is located in a much 
more diffuse H II region. NGC 6618PG is associated 
with both Ha em ission and ROSAT soft X-ray emis- 
sion ( Dunne et al.|[2003 ). The soft X-ray emission is also 
apparent in an archival Einstein observation taken 13 
years earlier than the i? 05^4 T survey. Significant differ- 
ences in morphology between these 2 images show that 
the X-ray emission is variable on ^10-yr timescales. We 

^® See |http:/ /www.mpifr-bonn.mpg.de/survey.htmll 
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Fig. 14.— MIPS 24 image highlighting the OB stars 
of NGC 6618PG (heavy diamonds) and the 90 cm radio ring 
(contours). Two candidate B stars associated with 24 /im 
emission are indicated by thin diamonds. Source numbers 
are from Table [S] 

therefore infer that the emission is dominated by unre- 
solved X-ray-bright point sources; these are the young 
stars of the NGC 6618PG cluster. High-resolution X-ray 
observations of the interior of M17 EB with Chandra or 
XMM-Newton are needed to separate individual cluster 
members from field stars and to detect any truly diffuse 
X-ray emission. 

The NGC 6618PG cluster must be significantly more 
evolved than NGC 6618 because, even in the absence of 
bright diffuse mid-IR emission from the interior of M17 
EB, we do not detect any candidate YSOs above the 
background level in the vicinity of the candidate ion- 
izing stars and X-ray emission. Given that elsewhere 
in the target field we detect at least a few solar-mass 
stars with optically thick disks, the absence of candi- 
date YSOs means that the age of NGC 6618PG is >2 
Myr (|Haisch. Lada. fc Lad"all2001al) . We therefore con- 
sider the possibility that the OB stars have evolved off 
the main sequence, and list possible spectral types of 
luminosity class HI for the most luminous stars, BD- 
16 4831 and BD-16 4826 (Nos. 17 and 18 in Table [5l). 
These stars place an upper limit on the age of NGC 
6618PG. Assuming negligible mass loss, an 06.5 HI star 
spent ~4.5 Myr as an 05 V star f ollowed by 5,0-5 Myr 
of post-main-sequence evolution (jBressan et al.l 119931 : 
iMartins. Schaerer. fc Hillieil 120051 ). for a maximum age 
of ^^5 Myr. Stars 17 and 18 are the best candidates 
for ionizing the H II region inside M17 EB, although we 
cannot rule out the possibility that other luminous stars 
once contributed to the ionization and have since gone 
supernova. 

Locations of all 5 candidate ionizing stars in NGC 
6618PG and the surrounding interior of M17 EB are 
marked on a MIPSGAL 24 /xm image in Figure [HI Also 
shown are co ntours of 90 cm radio continuum from VLA 
observations (jBrogan et al.l l2006l ) . An obvious feature 
of this image is a ring-shaped structure, centered at 
{l,b) = (15.23°, -0.72°) and -6' in diameter, that is 
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apparent in both the IR and radio continuum. While 
this structure is superficially suggestive of a SNR, careful 
analysis of available VLA data at 20 cm and 90 cm com- 
bined with 3.6 cm observations from the commissioning 
of the 100-m Robert C. Byrd Green Bank Telescope^^ 
reveals that the radio spectral index is consistent with 
thermal emission. There is no significant difference be- 
tween the spectral index of the radio ring and the rest 
of the M17 H II region. The 24 /xm morphology is remi- 
niscent of a stellar-win d bubble produced b y a luminous 
star (hke bubble N49 in lWatson et al.l|2008l ), but there is 
only a faint hint of a PDR in the GLIMPSE 8.0 ^m im- 
ages, and this structure is not detected in the HHT CO 
maps. We have fit SEDs to every Archive source located 
within the ring, and we find no candidate early-type star 
apart from those included in Figure fT4l and Table [3 Star 
17 is the closest O star to the center of the ring, but it 
is not central enough to explain the ring geometry. The 
origin of the ring is therefore puzzling. Recent structural 
models of the Ml 7 H II region suggest that the blister 
has broken to the North [Pellegrini e t al, 2007; PSC07), 
and hence could vent energy into the interior of M17 
EB. The ring structure may be material blown from the 
side of the M17 H II region by such venting. This raises 
the possibility that NGC 6618 contributes energy to M17 
EB. 

5.2. Constraints on the Expansion Timescale of M17 

EB 

M17 EB is a PDR around a faint, diffuse H II region. 
In the classical model for an expanding H II region (ne- 
glecting stellar winds), the timescale for expansion to a 
radius R is given by 
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from Equation 12-20 of lSpitzerl (|1978D . We assume typi- 
cal properties of ionized gas at = 8, 000 K: Isothermal 
sound speed Cu = 10 km s~^ and recombination co- 
efficient a^^) = 3.09 X 10~^^ cm^ s~^. The Stromgren 
radius Rs depends upon the ambient ISM gas density 
uq. We can make a rough estimate of no from the col- 
umn density of molecular gas swept up by the expan- 
sion of the bubble. Using the column density of molec- 
ular and atomic gas a.t v — 12-26 km s~^ derived from 
the CO emission (see §4.1), the total gas mass in M17 
EB (including M17 North and part of MC G15.9-0.7) is 
'-50,000 Mq. If we model M17 EB as a thin spherical 
shell with diameter 2R = 20 pc, then redistributing the 
gas throughout the shell interior yields a mean ambient 
density of uq 350 cm~'^, consistent with an inhomoge- 
neous molecular cloud environment. Photoionzed expan- 
sion driven by stars 17 and 18 in NGC 6618PG cannot 
produce an H II region the size of M17 EB within the 
stellar lifetimes (tphot > 5 Myr). 

The expansion of M17 EB could instead be dominated 
by stellar winds. A ccording to the analytical model of 
IWeaver et al.l (|1977t ). the timescale for expansion of a 
windblown bubble to radius R is 
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See |http: / /wiki. gb.nrao.edu/bin/view/Observing/M17publicImagi 



Typical wind luminosities for massive stars are 
iwind/iboi = 2.5 X 1 0~'^ for dwarfs and LwindZ -Lboi ~ 
5 X 10-3 for giants (|Vink. de Koder. k LameTi l2001h . 
For a given no, wind-driven expansion produces a large 
bubble mu ch more quic k ly th an photoionized expan- 
sion. The IWeaver etall ()1977[ ) model with no 350 
cm^-^ and the same 2 driving stars produces M17 EB in 
iwind < 1 Myr. 

The analytical models for both photoionized expansion 
and wind-driven expansion are oversimplifications of the 
physics governing the structure of H II regions. Both 
models underestimate the expansion timescales because 
they neglect radiative cooling by dust and, in the case of 
expansion driven by gas pressure (Equation 3), depres- 
surization caused by breaks in the bubble. Given the 
above caveats, the ideal case of photoionized expansion 
gives a lower limit on tphot that exceeds the stellar life- 
times. In the wind-driven case, unless radiative cooling 
increases twind by a factor >10, expansion driven by wind 
momentum (Equation 4) from the 2 most luminous stars 
in NGC 6618 PC can produce the bubble in iwind < 5 
Myr. In the other limiting case of twind > 2 Myr (see §5.1 
above), the average expansion velocity of the bubble is 
4.5 km s~^, and the current velocity is lower because 
the expansion slows with time. The signature of such 
slow expansion would be undetectable in our spatially- 
integrated CO line profiles (§4). 

5.3. Sequence of Events 

We now reconstruct a plausible star formation history 
for M17 based upon the wealth of observational data pre- 
sented in this work. 

Initially, the area shown within the white boundaries in 
Figure [15] was spanned by a large GMC complex that in- 
cluded the gas now observed in the M17 molecular cloud 
and MC G15.9-0.7. Between 2 and 5 Myr ago, a massive 
cloud subcomponent in the center of the complex under- 
went gravitational collapse, forming the first cluster of 
perhaps 2,000-3,000 stars, NGC 6618PG, the predeces- 
sor and possible progenitor of subsequent star formation 
in the complex. The OB stars in NGC 6618PG {blue di- 
amonds in Fig.ll5p ionized an H II region that expanded, 
sweeping ambient ISM gas and dust into a bubble, M17 
EB. 

The M17 molecular cloud, with initial mass >7 x 10* 
Mq, lay to the south of NGC 6618PG, in the path of 
the expanding M17 EB. What happened next is debat- 
able. The onset of star formation in M17 may have been 
triggered when M17 EB collided with the M17 molecu- 
lar cloud. But triggered star formation is impossible to 
prove, because we cannot prove that the molecular cloud 
would not have collapsed on its own, independent of the 
existence of NGC 6618PG. Indeed, it is curious that NGC 
6618, the more massive cluster, formed later than NGC 
6618PG. Measurements of the magnetic field strength 
in the PDR interface between the H II region and M17 
South (see Fig. [TS]) indicate that M17 South IS m ap- 
proximate dynamic equilibrium, with thermal, turbulent, 
and magnetic pressure combined supporting t he cloud 
(jBrogan fc Trolandll200ll: IPellegrini etal]|2007t ). These 
support mechanisms might have allowed M17 initially to 
resist gravitational collapse. 
The M17 molecular cloud began to collapse and form 
ie] stars not more than 2 Myr ago. We know that NGC 
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Fig. 15. — Overview of the extended massive young stellar population associated with M17 described in this work. Cyan 
contours of ^'^CO emission at u = 19 km s~^ are overlaid on a GLIMPSE 8 fim image. Candidate YSOs associated with the 
extended M17 molecular cloud and M17 EB are marked by colored circles as in Fig. 1101 Diamonds mark the known O stars 
of NGC 6618 (white), the 5 most massive OB stars detected in NGC 6618 (blue), and the 5 candidate OB stars identified by 
associated diffuse 24 fim emission (orange). 



6618 is younger than NGC 6618PG, because: (1) Despite 
being a more luminous cluster exhibiting many signs 
of strong stellar winds (|Townslev et all 120031: IPSCOTI : 
iPovich et a l. 2008), NGC 6618, unlike NGC 6618PG, has 
not yet dispersed its natal cloud; and (2) NGC 6618 
has numerous highly-embedded intermediate- to high- 
mass stars (ongoing massive star formation), while NGC 
6618PG has no intermediate-mass stars with optically 
thick disks (massive star formation has ceased). The 
burst of star formation that created the most massive 
members of NGC 6618 was probably very rapid. The 
time-averaged SFR in the central regions of M17 has been 
0.008-0.01 Mq yr-i. 

The presence of NGC 6618PG explains some puzzling 
features of the classic "blister" structure of the M17 
H II region. Why should a blister H II region form in 
the first place? This question has received its share of 
theoretical attention in recent years. Two-dimensional 
models of molecular clouds with uniform density show 
that gravitational acceleration p eaks at the edges of 
clouds, sweeping material in ward (|Burkert fc HartmannI 
I2004f) . iHeitsch et all ()2008D present 3-dimensional nu- 
merical simulations that lead to the formation of mas- 
sive protostellar cores on the edges of molecular clouds 
facilitated by gravitational instabilities that focus large- 
scale colliding flows. M17 could be a good example of 
the end-product of such a process, with its cluster of 
very massive young stars forming on the northern side 
of the massive M17 South cloud. The CO maps show a 
narrow molecular ridge tracing M17 North that appears 



to be pressure-bounded on 2 sides. The ridge hes be- 
tween NGC 6618 and NGC 6618PG (contours between 
NGC 6618 and NGC 6618PG in Fig. [IS]). The source 
of outward pressure from the stars ionizing the M17 H II 
region is well-known, and we have now identified a poten- 
tial source of counteracting pressure from the expansion 
of M17 EB. This provides a natural explanation for the 
morphology of the M17 North cloud and the associated 
northern half of the M17 H II region and PDR. 

As the Ml 7 H II region pushes outward from M17 
South, M17 North is eroded. Once this barrier is 
breached, the built-up pressure in the H II region can be 
released into the diffuse interior of M17 EB. There are 
observational clues that this may alread y have occurred. 
The d iffuse X-ray emission identified by iTownslev et al.l 
(|2003f ) traces hot plasma filling the interior cavity of the 
M17 H II region. This plasma appears to spill out of 
the H II region to the East and North (Townsley et al. 
in prep.). The thermal radio/IR ring shown in Figure 
[14] may also be evidence of a blowout from the northern 
side of the H II region. In this scenario, it is possible that 
NGC 6618 now pumps energy into M17 EB, even though 
the bubble was originally driven by NGC 6618PG. Addi- 
tional observations in X-rays and optical recombination 
lines are needed to reveal the detailed structure of Ml 7 
EB. 

The expansion of M17 EB is asymmetric, with NGC 
6618PG offset from the center of the bubble toward its 
southern rim, near M17 itself. This is reasonable given 
that the presence of the large molecular mass in M17 in- 
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dicates a strong density gradient across the bubble. Any 
energy input from NGC 6618 into M17 EB would provide 
an additional push to expand the bubble preferentially 
away from the M17 H II region. 

The northward expansion of Ml 7 EB overtook another 
molecular cloud, MC 15.9-0.7. The central regions of MC 
G15.9-0.7 lack IR signatures of star formation. As Figure 
[TSl shows, only the side of this cloud that intersects the 
8 /im emission from M17 EB exhibits a significant con- 
centration of candidate YSOs after the removal of back- 
ground contaminants. Hence the circumstantial evidence 
for triggered star formation is stronger here than in the 
case of NGC 6618. Outside of M17 itself, the bubble- 
cloud interface produces the brightest CO (J = 2 ^ 1) 
emission at w = 20 km s~^. The current SFR in the 
molecular gas swept up by M17 EB and in the extended 
outer regions of M17 is >16% of the SFR that produced 
the well-studied young stellar population of the M17 H II 
region. 

6. CONCLUSIONS 

We have used IR, radio, and X-ray survey data com- 
bined with targeted millimeter observations to show that 
star formation in the M17 complex extends ~0.5° far- 
ther to the north than previously thought. Our ^^CO 
and ^■^CO line maps cover a 0.72 deg^ area that contains 
1.35 X 10^ Mq in gas aX v = 12-26 km 8"^ The mor- 
phology of the CO map (Fig. is dominated by an 
extended bubble, M17 EB, occupying the space between 
the well-studied M17 molecular cloud and a neighboring 
molecular cloud to the north, MC G15.9-0.7. 

We have found 406 candidate YSOs in a 1.5° x 1° tar- 
get field centered on M17 EB by fitting model SEDs to 
fluxes of sources in the GLIMPSE Archive supplemented 
by MIPSGAL and MSX photometry. We characterized 
the contamination from unassociated sources in our YSO 
sample by analyzing a control field exhibiting no obvious 
signs of massive star formation. We were able to identify 
and remove the brightest >50% of dust-rich AGB stars 
from our sample. We conclude that >80% of sources in 
the control field sample are candidate YSOs distributed 
at unknown distances along the long sightline through 
the inner Galactic plane. 

We created a map of the extinction produced by the 
molecular cloud structures associated with M17. Out- 
side of a few small regions of high gas column density, 
across ~90% of our target field the extinction is too low 
to significantly reduce mid-IR detections of background 
sources. The complex, spatially varying extinction in the 
M17 molecular cloud itself creates special challenges for 
observations of the young stellar population associated 
with NGC 6618. The average extinction through the 
cloud cavity occupied by the central M17 H II region is 
Av < 17 mag, too low to prevent the detection of a sig- 
nificant fraction of background stars even for relatively 
shallow near-IR observations. Recent results reported by 
H08, in particular the high IR excess fraction and young 
age of 0.5 Myr for the young stellar population in M17, 
appear to have been compromised by the use of an IR- 
selected sample that was not corrected for background 
contamination. 

Approximately half of the candidate YSOs in our tar- 
get field are unassociated contaminants. In studies of 
massive young clusters, IR excess sources are commonly 



assumed to be associated by virtue of their apparent 
youth, but this assumption is not valid in the case of star- 
formation regions near the inner Galactic mid-plane. By 
selecting candidate YSOs exhibiting a significant degree 
of clustering with respect to the control sample, we suc- 
cessfully removed the predicted level of contamination. 

Among 195 candidate YSOs overlapping the area oc- 
cupied by the extended M17 complex, 96 have high sta- 
tistical likelihood of association with M17 once the con- 
tribution from unassociated foreground and background 
sources is removed. Modeling the YSO population, we 
find that our sample consists pri marily of intermediate- 
mass YSOs with disks (Stage II; iRobitaille erani2006l ) 
and highly embedded (Stage O/I) YSOs. We model the 
accretion rates for the Stage O/I sources and find that 
their median age is <0.075 Myr, leading us to derive a 
circumstellar disk lifetime of ~0.5 Myr for sources with 
> 3 M0 (corresponding to B-type stars on the main 
sequence). Candidate YSOs distributed around the rim 
of M17 EB and the extended outer regions of M17 rep- 
resent a >16% addition to the star formation rate that 
has produced the young stellar population in the central 
regions of M17, as observed by Chandra (BFT07). The 
concentration of star formation activity tracing the in- 
terface between M17 EB and MC G15.9-0.7 appears to 
have been triggered by the expansion of the bubble. 

We have identified 5 candidate ionizing stars inside 
M17 EB. By fitting stellar atmosphere models to their IR 
SEDs we constrained the luminosity distance of these 5 
sources, showing that their IR luminosities give OB spec- 
tral types at the M17 distance. An associated lower-mass 
population of young stars is suggested by X-ray survey 
observations by the Einstein and i? 05*^ T" satellites, lead- 
ing us to deduce the presence of a young cluster, NGC 
6618PG. The two earhest-type stars in NGC 6618PG, 
BD-16 4831 and BD-16 4826 (stars 17 and 18 in Table 
[5|), have luminosities equivalent to main-sequence spec- 
tral types 04-05, but the stars may be evolved. The 
stellar winds from these two stars carry sufficient mo- 
mentum to produce an H II region the size of M17 EB 
in <5 Myr. The current expansion velocity of M17 EB 
is low, a few km s~^ maximum. 

We also fit stellar atmosphere SEDs to 11 of the prin- 
cipal ionizing stars in NGC 6618. At the newly-revised 
2.1 kpc distance (H08), we find that a high fraction of 
the O stars in NGC 6618 are too luminous in the mid-IR 
to be single stars on the main sequence, yet the stars are 
too young to have evolved off the main sequence. This 
finding supports the hypothesis of H08 that most of the 
O stars in M17 are in binary systems with approximately 
equal-mass components. 

The wealth of observational data paints a picture of 
propagating star formation in this region. Three main 
waves of star formation can be identified: (1) The for- 
mation of the NGC 6618PG cluster (2,000-3,000 stars), 
2-5 Myr in the past, followed by the expansion of an 
H II region delineated by M17 EB; (2) The rapid col- 
lapse of the M17 molecular cloud (current mass >6 x 10** 
Mq) within the last 2 Myr, producing the massive NGC 
6618 cluster and associated massive star formation that 
continues to the present day (11,000-14,000 stars); and 
(3) The onset of star formation in MC G15.9-0.7 and 
the extended outer regions of M17 (>1,000 stars). The 
possibility that the later waves of star formation were 
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triggered by the first raises a tantalizing question: Can 
a massive progenitor cluster trigger the formation of an 
even more massive daughter cluster? 

Because they are beacons observable at large dis- 
tances across the local universe, H II regions are com- 
monly used as SFR tracers in external galaxies, and have 
been used to estimate the SFR of the Milky Way as 
well ([Smith. Biermann. fc Mezgeilll978l ). The canonical 
Galactic SFR of a few M0 yr~^ is equivalent to sev- 
eral hundred H II regions the size of M17. Observa- 
tional tracers tied to H II regions are sensitive only to 
the most massive <1% of stars, and the uncertainty in- 
volved in extrapolating over the IMF to the low-mass 
stars that make up the bulk of the stellar mass is enor- 
mous. Star formation distributed at lower densities but 
over large spatial volumes is increasingly observed in 
promin ent Galactic molecula r cloud complexes, includin g 
Orion (|Megeath et al.l[200l . W5 (|Koemg et al.l[200l . 
and now M17. It is not yet clear if low-density star for- 
mation outside of bright H II regions or dense GMCs 
significantly increases the galactic SFRs derived from ob- 
servational tracers of ionized or molecular gas. 
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